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Abstract 

We study the phase diagram of Q-state Potts models, for Q = 4cos 2 (-7r/p) a Beraha 
number (p > 2 integer), in the complex-temperature plane. The models are defined 
onlxJV strips of the square or triangular lattice, with boundary conditions on the 
Potts spins that are periodic in the longitudinal (./V) direction and free or fixed in the 
transverse (L) direction. The relevant partition functions can then be computed as 
sums over partition functions of an A p _i type RSOS model, thus making contact with 
the theory of quantum groups. We compute the accumulation sets, as N — > oo, of 
partition function zeros for p = 4, 5, 6, oo and L = 2, 3, 4 and study selected features 
for p > 6 and/or L > 4. This information enables us to formulate several conjectures 
about the thermodynamic limit, L — ► oo, of these accumulation sets. The resulting 
phase diagrams are quite different from those of the generic case (irrational p). For 



free transverse boundary conditions, the partition function zeros are found to be dense 
in large parts of the complex plane, even for the Ising model (p = 4). We show how 
this feature is modified by taking fixed transverse boundary conditions. 

Key Words: Potts model, RSOS model, Beraha number, limiting curve, quantum groups 
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1 Introduction 



The Q-state Potts model [1,2] can be denned for general Q by using the Fortuin-Kasteleyn 
(FK) representation [3,4]. The partition function Z G (Q;v) is a polynomial in the variables 
Q and v. This latter variable is related to the Potts model coupling constant J as 



v = e J -l. (1.1) 



It turns out useful to define the temperature parameter x as 



x = 7Q (L2) 



and to parameterize the interval Q G (0, 4] as 



Q = 4 cos y-J , p G (2,oo]. (1.3) 

For generic 1 values of Q, the main features of the phase diagram of the Potts model 
in the real (Q,x)-plane have been known for many years [2,5]. It contains in particular a 
curve xfm{Q) > of ferromagnetic phase transitions which are second-order in the range 
< Q < 4, the thermal operator being relevant. The analytic continuation of the curve 
xfm{Q) into the antiferromagnetic regime yields a second critical curve xbk(Q) < with 
< Q < 4 along which the thermal operator is irrelevant. Therefore, for a fixed value of Q, 
the critical point xbk(Q) acts as the renormalization group (RG) attractor of a finite range 
of x values: this is the Berker-Kadanoff (BK) phase [6,7]. 

The generic phase diagram is shown in Fig. 1. Since the infinite-temperature limit (x = 0) 
and the zero-temperature ferromagnet (\x\ = oo) are of course RG attractive, consistency of 
the phase diagram requires that the BK phase be separated from these by a pair of RG repul- 
sive curves x±(Q) < 0. The curve x+(Q) is expected to correspond to the antiferromagnetic 
( AF) phase transition of the model [8] . 

The above scenario thus essentially relies on the RG attractive nature of the curve x B k(Q), 
and since this can be derived [5] from very general Coulomb gas considerations, the whole 
picture should hold for any two-dimensional lattice. But it remains of course of great interest 
to compute the exact functional forms of the curves xfm(Q), xbk(Q), an d x±(Q) — and the 
corresponding free energies — for specific lattices. 

The square-lattice Potts model is the best understood case. Here, Baxter [2,9] has found 
the exact free energy along several curves x = x c (Q): 

+ 1 (FM) 

~ J k + \lW ( AF ) 
Xc (Q) = ) ^ V Q (1 .4) 

* -1 (BK) 

2 U-Q 

To ~ ■ — 



(AF) 



where x c — 1 and x c = — 1 can be identified respectively with xfm(Q) and xbk(Q)- The 
curves x± = —2/y/Q ± -^/(4 — Q)/Q are mutually dual (and hence equivalent) curves of AF 



More precisely, a "generic" value of Q corresponds to an irrational value of the parameter p defined in 
Eq. (1.3). This point will be made more precise in Section 2 below. 
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phase transitions, which are again second-order in the range < Q < 4. These curves also 
form the boundaries of the x- values controlled by the BK fixed point [7] , as outlined above. 
Note that the four points x c (q) in Eq. (1.4) correspond to the points where the circles 

\x\ = 1 (1.5a) 
-4 = kI^- (1.5b) 



Q 

cross the real x-axis. These two circles intersect at the points 

x = - e ±t7T/p (1.6) 

which will be shown below to play a particular role in the phase diagram (see Conjec- 
ture 4.1.1). 

In the case of a triangular lattice, Baxter and collaborators [10-12] have found the free 
energy of the Potts model along the curves 

^/Qx 3 + 3x 2 = 1, (1.7a) 

* = -ja- (17b) 



The upper branch of Eq. (1.7a) is identified with the ferromagnetic critical curve £fm(<3)- 
We have numerical evidence that the middle and lower branches correspond respectively 
to xbk(Q) and x-{Q), the lower boundary of the BK phase. The position of x+{Q), the 
upper branch of the BK phase, is at present unknown [13] (but see Ref. [14] for the Q — > 
limit). Along the line (1.7b) the Potts model reduces to a coloring problem, and the partition 
function is here known as the chromatic polynomial. The line (1.7b) belongs to the RG basin 
of the BK phase for < Q < 2 + ^ [15]. 

The critical properties — still with Q taking generic values — for these two lattices are to 
a large extent universal. This is not so surprising, since the critical exponents can largely be 
obtained by Coulomb gas techniques (although the antiferromagnetic transition still reserves 
some challenges [8]). Thus, there is numerical evidence that the exponents along the curves 
Xfm(Q), %bk(Q) and X-(Q) coincide, whereas the evidence for the curve x + (Q) is non- 
conclusive [14]. On the other hand, on the less-studied triangular lattice we cannot yet 
exclude the possible existence of other curves of second-order phase transitions that have no 
counterpart on the square lattice. 

But in general we can only expect universality to hold when the Boltzmann weights in 
the FK representation are non-negative (i.e., for Q > 0, v > 0), or when the parameter p 
takes generic (i.e., irrational) values. The present paper aims at studying the situation when 
p takes non-generic values; for simplicity we limit ourselves to the case of integer p > 2. The 
number of spin states is then equal to a so-called Beraha number B p 



Q = B p = 4cos 2 (Jj , p = 3,4,5,.... (1.8) 

The special physics at rational values of p is intimately linked to the representation theory 
of the quantum group C/ ? (577(2)), the commutant of the Temperley-Lieb algebra, when the 
deformation parameter q is a root of unity. As we shall review in Section 2 below, the 
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quantum group symmetry of the Potts model at rational p implies that many eigenvalues of 
the transfer matrix in the FK representation have zero amplitude or cancel in pairs because 
of opposite amplitudes; these eigenvalues therefore become spurious and do not contribute 
to the partition function [6,7]. 

Remarkably, for p integer and x inside the BK phase, even the leading eigenvalue acquires 
zero amplitude. Moreover, all the eigenvalues which scale like the leading one in the thermo- 
dynamic limit vanish from the partition function, and so, even the bulk free energy f(p;x) 
is modified [8]. In other words, f(p;x) experiences a singularity whenever p passes through 
an integer value. This means in particular that for p integer the critical behavior can either 
disappear, or be modified, or new critical points (and other non-critical fixed points) can 
emerge. 

For the sake of clarity, we discuss the simplest example of this phenomenon. Consider, 
on the square lattice, on one hand the Q — > 2 state model (i.e., with Q tending to 2 through 
irrational values of p) and on the other the Q = 2 Ising model (i.e., with fixed integer p — 4). 
For the former case, the generic phase diagram and the associated RG flows are shown in 
the top part of Fig. 2. The three critical points xfm and x± have central charge c = 1/2, 
while the fourth one xbk has c = —25/2. For the latter case, new non-critical fixed points 
appear (by applying the duality and Z 2 gauge symmetries to the one at x — 0), and the 
RG flows become as shown in the bottom part of Fig. 2. One now has c = 1/2 for all four 
critical fixed points. (We shall treat the Ising model in more detail in Section 7.1 below.) 

By contrast to the universality brought out for generic Q, the phase diagram and critical 
behavior for integer p is likely to have lattice dependent features. Let us give a couple of 
examples of this non-universality. The zero-temperature triangular-lattice Ising antiferro- 
magnet, (Q,v) = (2,-1), is critical and becomes in the scaling limit a free Gaussian field 
with central charge c = 1 [16-18], whereas the corresponding square-lattice model is non- 
critical, its partition function being trivially Z = 2. While this observation does not in 
itself imply non-universality, since the critical temperature is expected to be lattice depen- 
dent (as is the value of £fm(Q)), the point to be noticed is that for no value of v does the 
Q = 2 square-lattice model exhibit c = 1 critical behavior. In the same vein, the square- 
lattice Potts model with (Q,v) = (3,-1) is equivalent to a critical six-vertex model (at 
A = 1/2) [19,20], with again c = 1 in the scaling limit, whereas now the corresponding 
triangular-lattice model is trivial (Z = 3). Now, the triangular- lattice model does in fact 
exhibit c = 1 behavior elsewhere (for x = but the compactification radius is different 
from that of the square-lattice theory and accordingly the critical exponents differ. Finally, 
(Q,v) = (4, —1) is a critical c = 2 theory on the triangular lattice [21,22], but is non-critical 
on the square lattice [23] . 

Because of the eigenvalue cancellation scenario sketched above, the FK representation 
is not well suited 2 for studying the Potts model at integer p. Fortunately, for Q = B p 
there exists another representation of the Potts model, in terms of an RSOS model of the 
Ap-i type [24], in which the cancellation phenomenon is explicitly built-in, in the sense that 
for generic values of x all the RSOS eigenvalues contribute to the partition function. On 
the square lattice, the RSOS model has been studied in great detail [24-27] at the point 

2 We here tacitly assume that the study relies on a transfer matrix formulation. This is indeed so in most 
approaches that we know of, whether they be analytical or numerical. An exception would be numerical 
simulations of the Monte Carlo type, but in the most interesting parts of the phase diagram this approach 
would probably not be possible anyway, due to the presence of negative Boltzmann weights. 
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x = xfm = 1, where the model happens to be homogeneous. Only very recently has the case 
of general real x ^ 1 (where the RSOS model is staggered, i.e., its Boltzmann weights are 
sublattice dependent) attracted some attention [8], and no previous investigation of other 
lattices (such as the triangular lattice included in the present study) appears to exist. 

The very existence of the RSOS representation has profound links [27,28] to the represen- 
tation theory of the quantum group U q (SU(2)) where the deformation parameter q defined 
by 

Q = (q + q- 1 ) 2 = B p , q = exp(ivr/p) , (1.9) 

is a root of unity. To ensure the quantum group invariance one needs to impose periodic 
boundary conditions along the transfer direction. Further, to ensure the exact equivalence 
between Potts and RSOS model partition functions the transverse boundary conditions must 
be non-periodic. 3 For definiteness we shall therefore study square- or triangular-lattice strips 
of size L x N spins, with periodic boundary conditions in the iV-direction. The boundary 
conditions in the L-direction are initially taken as free, but we shall later consider fixed 
transverse boundary conditions as well. For simplicity we shall henceforth refer to these 
boundary conditions as free cyclic and fixed cyclic. 4 

Using the RSOS representation we here study the phase diagram of the Potts model at 
Q = B p through the loci of partition function zeros in the complex x-plane. According to 
the Beraha-Kahane- Weiss theorem [30], when N — > oo, the accumulation points of these 
zeros form either isolated limiting points (when the amplitude of the dominant eigenvalue 
vanishes) or continuous limiting curves Bl (when two or more dominant eigenvalues become 
equimodular); we refer to Ref. [32] for further details. In the RSOS representation only 
the latter scenario is possible, since all amplitudes are strictly positive. 5 As usual in such 
studies, branches of Bl that traverse the real x-axis for finite L, or "pinch" it asymptotically 
in the thermodynamic limit L — > oo, signal the existence of a phase transition. Moreover, 
the finite-size effects and the impact angles [33] give information about the nature of the 
transition. 

The limiting curves B L constitute the boundaries between the different phases of the 
model. Moreover, in the present set-up, each phase can be characterized topologically by 
the value of the conserved quantum group spin S z , whose precise definition will be recalled 
in Section 2 below. (A similar characterization of phases of the chromatic polynomial was 
recently exploited in Ref. [34], but in the FK representation). One may think of S z as a kind 
of "quantum" order parameter. A naive entropic reasoning would seem to imply that for 
any real x the ground state (free energy) has S z = 0, since the corresponding sector of the 
transfer matrix has the largest dimension. It is a most remarkable fact that large portions 
of the phase diagram turn out have S z ^ 0. 

We have computed the limiting curves Bl in the complex x-plane completely for p = 

3 There are however some intriguing relationships between modified partition functions with fully periodic 
boundary conditions [29]. We believe that the RSOS model with such boundary conditions merits a study 
similar to the one presented here, independently of its relation to the Potts model. 

4 It is convenient to introduce the notation Lp x Np (rcsp. Lx x Np) for a strip of size L x N spins with 
free (resp. fixed) cyclic boundary conditions. 

5 Sokal [31, Section 3] has given a slight generalization of the Beraha-Kahane- Weiss theorem. In partic- 
ular, when there are two or more equimodular dominant eigenvalues, the set of accumulation points of the 
partition-function zeros may include isolated limiting points when all the eigenvalues vanish simultaneously. 
See Section 3.1.1 for an example of this possibility. 
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4, 5,6, oo and L = 2,3,4 for both lattices. Moreover, we have studied selected features 
thereof for p > 6 and/or L > 4. This enables us to formulate several conjectures about the 
topology of Bl which are presumably valid for any L, and therefore, provides information 
about the thermodynamic limit L — > oo. The resulting knowledge is a starting point for 
gaining a better understanding of the fixed point structure and renormalization group flows 
in these Potts models. 

Our work has been motivated in particular by the following open issues: 

1. As outlined above, the eigenvalue cancellation phenomenon arising from the quantum 
group symmetry at integer p modifies the bulk free energy in the Berker-Kadanoff 
phase. For the Ising model we have seen that this changes the RG nature (from 
attractive to repulsive) of the point well as its critical exponents (from c = 
—25/2 to c = 1/2). But for general integer p it is not clear whether x^k will remain a 
phase transition point, and assuming this to be the case what would be its properties. 

2. The chromatic line x = — \j\[Q does not appear to play any particular role in the 
generic phase diagram of the square-lattice model. By contrast, it is an integrable 
line [11,12] for the generic triangular-lattice model. Qua its role as the zero-temperature 
antiferromagnet one could however expect the chromatic line to lead to particular (and 
possibly critical) behavior in the RSOS model. Even when critical behavior exists in 
the generic case (e.g., on the triangular lattice) the nature of the transition may change 
when going to the case of integer p (e.g., from c = —25/2 for the Q — > 2 model to c = 1 
for the zero-temperature Ising antiferromagnet). 

3. Some features in the antiferromagnetic region might possibly exhibit an extreme depen- 
dence on the boundary conditions, in line with what is known, e.g., for the six- vertex 
model. It is thus of interest to study both free and fixed boundary conditions. To give 
but one example of what may be expected, we have discovered — rather surprisingly— 
that with free cyclic boundary conditions the partition function zeros are actually dense 
in substantial parts of the complex plane: this is true even for the simplest case of the 
square-lattice Ising model. 

4. A recent numerical study [8] of the effective central charge of the RSOS model with 
periodic boundary conditions, as a function of x, has revealed the presence of new 
critical points inside the BK phase. In particular, strong evidence was given for a 
physical realization of the integrable flow [35] from parafermion to minimal models. 
The question arises what would be the location of these new points in the phase 
diagram. 

5. In the generic case, the spin S z of the ground state may be driven to arbitrary large 
values upon approching the point (Q,x) = (4,-1) from within the BK phase [7,34]. 
Is a similar mechanism at play for integer p? 

The paper is organized as follows. In Section 2 we introduce the RSOS models and 
describe their precise relationship to the Potts model, largely following Refs. [24,27,28]. We 
then present, in Section 3, the limiting curves found for the square-lattice model with free 
cyclic boundary conditions, leading to the formulation of several conjectures in Section 4. 
Sections 5-6 repeat this programme for the triangular-lattice model. In Section 7 we discuss 
the results for free cyclic boundary conditions, with special emphasis on the thermodynamic 
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limit, and motivate the need to study also fixed cyclic boundary conditions. This is then 
done in Sections 8-9. Finally, Section 10 is devoted to our conclusions. An appendix gives 
some technical details on the dimensions of the transfer matrices used. 

2 RSOS representation of the Potts model 

The partition function of the two-dimensional Potts model can be written in several 
equivalent ways, though sometimes with different domains of validity of the relevant param- 
eters (notably Q). The interplay between these different representations is at the heart of 
the phenomena we wish to study. 

The spin representation for Q integer is well-known. Its low-temperature expansion gives 
the FK representation [3, 4] discussed in the Introduction, where Q is now an arbitrary 
complex number. The (interior and exterior) boundaries of the FK clusters, which live on 
the medial lattice, yield the equivalent loop representation with weight Q 1 ^ 2 per loop. 

An oriented loop representation is obtained by independently assigning an orientation to 
each loop, with weight q (resp. q^ 1 ) for counterclockwise (resp. clockwise) loops, cf. Eq. (1.9). 
In this representation one can define the spin S z along the transfer direction (with paral- 
lel/antiparallel loops contributing ±1/2) which acts as a conserved quantum number. Note 
that S z = j means that there are at least j non-contractible loops, i.e., loops that wind around 
the periodic (N) direction of the lattice. The weights q ±l can be further redistributed locally, 
as a factor of q a / 2n for a counterclockwise turn through an angle a [2]. While this redistri- 
bution correctly weights contractible loops, the non-contractible loops are given weight 2, 
but this can be corrected by twisting the model, i.e., by inserting the operator q Sz into the 
trace that defines the partition function. 

A partial resummation over the oriented-loop splittings at vertices which are compatible 
with a given orientation of the edges incident to that vertex now gives a six-vertex model 
representation [36]. Each edge of the medial lattice then carries an arrow, and these arrows 
are conserved at the vertices: the net arrow flux defines S z as before. The six-vertex model 
again needs twisting by the operator q Sz to ensure the correct weighing in the S z ^ sectors. 
The Hamiltonian of the corresponding spin chain can be extracted by taking the anisotropic 
limit, and is useful for studying the model with the Bethe Ansatz technique [2]. The fact 
that this Hamiltonian commutes with the generators of the quantum group U q (SU(2)) links 
up with the nice results of Saleur and coworkers [6,7,27,28]. 

Finally, the RSOS representation [24, 27, 28] emerges from a certain simplification of the 
above representations when q = exp(in/p) is a root of unity (see below). 

All these formulations of the Potts model can be conveniently studied through the cor- 
responding transfer matrix spectra: these give access to the limiting curves Bl, correlation 
functions, critical exponents, etc. 

In the FK representation the transfer matrix Tp^(L) is written in a basis of connectivities 
(set partitions) between two time slices of the lattice (see Ref. [34] for details), and the 
transfer matrix propagates just one of the time slices. Each independent connection between 
the two slices is called a bridge; the number of bridges j is a semi-conserved quantum number 
in the sense that it cannot increase upon action of the transfer matrix. The bridges serve to 
correctly weight the clusters that are non-contractible with respect to the cyclic boundary 



8 



conditions. 6 This is accomplished by writing the partition function as 

Z FK = (f\T^(Lf\ i ) = J2^^ (2.1) 

i>i 

for suitable initial and final vectors \i) and (/|. The vector \i) identifies the two time slices, 
while (f\ imposes the periodic boundary conditions (it "reglues" the time slices) and weighs 
the resulting non-contractible clusters. Note that these vectors conspire to multiply the 
contribution of each eigenvalue Aj by an amplitude ctj = (Xi(Q): this amplitude may vanish 
for certain values of Q. 

On the other hand, in the six-vertex representation the transfer matrix is written in the 
purely local basis of arrows, whence the partition function can be obtained as a trace (which 
however has to be twisted by inserting q Sz as described above). But even without the twist 
the eigenvalues are still associated with non-trivial amplitudes, as we now review. 

Let us consider first a generic value of q, i.e., an irrational value of p. The U q (SU(2)) 
symmetry of the spin chain Hamiltonian implies that one can classify eigenvalues according 
to their value j of S z , and consider only highest weights of spin j. Define now K 1>2 j+i(p, L; x) 
as the generating function of the highest weights of spin j, for given values of p, L and x. The 
partition function of the untwisted six-vertex model with the spin S (not S z ) fixed to j is 
therefore (2j + 1) Ki 2j+i(p, L; x). Imposing the twist, the corresponding contribution to the 
partition function of the Potts model becomes Sj(p) K^j+iip, L; x), where the g-deformed 
number Sj(p) = (2j + l) q is defined as follows 

sm(7r/p) 

S has a simple interpretation in the FK representation as the number of bridges, whereas it is 
S z which has a simple interpretation in the six- vertex model representation as the conserved 
current. 

Different representations correspond to choosing different basis states: a given cluster 
state is an eigenvector of S, but not S z , and a given vertex state is an eigenvector of S z , but 
not S. The eigenvectors of the Hamiltonian are eigenvectors of both S and S z , and are thus 
combinations of vertex states (or of cluster states if one works in the FK representation). 
But note that the dimensions of the transfer matrix are not exactly the same in the vertex 
and the FK representations, as the 2j + 1 possible values of S z for a given S — j are not taken 
into account in the same way: in the vertex representation, it corresponds to a degeneracy 
of the eigenvalues, whereas in the FK representation it appears because of the initial and 
final vectors which sandwich the transfer matrix in Eq. (2.1). 

The total partition function of the Q-state Potts model on a strip of size Lp x Np can 
therefore be exactly written as [24, 27, 28] 

L 

(Q; v) = Q LN ' 2 Sj(p) K 1>2j+1 (p, L; x) (2.3) 

3=0 

Note that the summation is for < j < L, as the maximum number of bridges is equal to 
the strip width L. 

6 In particular, the restriction of Tp^(L) to the zero- bridge sector is just the usual transfer matrix Tfk in 
the FK representation, i.e., the matrix used in Ref. [32] to study the case of fully free boundary conditions. 
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For p rational, Eq. (2.3) is still correct, but can be considerably simplified. In the context 
of this paper we only consider the simplest case of p integer. Indeed, note that using Eq. (2.2), 
we obtain that, for any integer n, 

S {n+ i )p -i-j{p) = -Sj(p) (2.4a) 
S np+j (p) = Sj(p). (2.4b) 

Therefore, after factorization, Eq. (2.3) can be rewritten as 

L(P"2)/2J 

{ q. v) = qLN/2 J- S,(p)xi,2 H i(p,L;x), (2.5) 

3=0 

where 

Xi,2j+i(p,L;x) = ^[K li 2( n p+j )+ i(p,L;x) - K 1: 2(( n+ i)p~i-j)+i(p,L;x)] . (2.6) 

n>0 

For convenience in writing Eq. (2.6) we have defined Ki t2 j+i(p, L; x) = for j > L. Note 
that the summation in Eq. (2.5) is now for < j < [(p — 2) / 2j . Furthermore, Xi,2j+i(p, L; x) 
is a lot simpler that it seems. Indeed, when p is integer, the representations of U q (SU (2)) mix 
different values of j related precisely by the transformations j — * j+np and j — > (n+l)p— 1— j 
[cf. Eq. (2.4)]. Therefore, a lot of eigenvalues cancel each other in Eq. (2.6). This is exactly 
why the transfer matrix in the FK representation contains spurious eigenvalues, and is not 
adapted to the case of p integer. 

The representation adapted to the case of p integer is the so-called RSOS representation. 
It can be proved that Xi,2j+i is the partition function of an RSOS model of the type [24] 
with given boundary conditions [27] (see below). In this model, heights hi — 1,2, ... ,p — 1 
are defined on the union of vertices and dual vertices of the original Potts spin lattice. 
Neighboring heights are restricted to differ by ±1 (whence the name RSOS = restricted 
solid-on-solid). The boundary conditions on the heights are still periodic in the longitudinal 
direction, but fixed in the transverse direction. More precisely, the cyclic strip Lp x Np has 
precisely two exterior dual vertices, whose heights are fixed to 1 and 2j + 1 respectively. It 
is convenient to draw the lattice of heights as in Figures 3-4 (showing respectively a square 
and a triangular-lattice strip of width L = 2), i.e., with N exterior vertices above the upper 
rim, and N exterior vertices below the lower rim of the strip: all these exterior vertices close 
to a given rim are then meant to be identified. 

For a given lattice of spins, the weights of the RSOS model are most easily defined by 
building up the height lattice face by face, using a transfer matrix. The transfer matrix 
adding one face at position i is denoted Hi = xli + (resp. V* = U + xti) if it propagates a 
height hi — > h\ standing on a direct (resp. a dual) vertex, where Jj = 5(hi, h'j) is the identity 
operator, and ej is the Temperley-Lieb generator in the RSOS representation [24]: 

c , =j , V] , M N!yf!^. (2 . 7 ) 

sm(irhj-i/p) 

Note that all the amplitudes Sj(p) entering in Eq. (2.5) are strictly positive. Therefore, 
for a generic value of the temperature x, all the eigenvalues associated with Xi,2j+i(p, L; x) 
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for < 2j + 1 < p contribute to the partition function. 7 This is the very reason why we 
use the RSOS representation. Recall that there are analogous results in conformal field 
theory [37]. In fact, for x equal to Xfm(Q) and in the continuum limit, K 12 j+i corresponds 
to the generating function of a generic representation of the conformal symmetry with Kac- 
table indices r = 1 and s = 2j + 1, whereas Xi,2j+i corresponds to the generating function 
(character) of a minimal model. Thus, Eq. (2.6) corresponds to the Rocha-Caridi equation 
[38], which consists of taking into account the null states. One could say that the FK 
representation does not identify all the states differing by null states, whereas the RSOS 
representation does. Therefore, the dimension of the transfer matrix is smaller in the RSOS 
representation than in the FK representation. 

The computation of the partition functions Xi,2j+i(p, L; x) can be done in terms of transfer 
matrices T 12?+1 , denoted in the following simply by T 2 j+i- In particular, for a strip of size 
L x N, we have that 

Xi,2j+i(p, L; x) = trJ 2j+1 (p,L;x) N (2.8) 

Note that this is a completely standard untwisted trace. The transfer matrix T 2j - + i(L; x) acts 
on the space spanned by the vectors \ho, h±, . . . , h 2 i), where the boundary heights ho — 1 
and hih = 2 j + 1 are fixed. The dimensionality of this space is discussed in Appendix A. 
For any fixed h and h 2 L, this dimensionality grows asymptotically like ~ Q L . 

Remarks. 1) Our numerical work is based on an automatized construction of T 2j+1 . To 
validate our computer algorithm, we have verified that Eq. (2.5) is indeed satisfied. More 
precisely, given Q = B p , and for fixed L and N, we have verified that 

Zl f *n p (Q;v) = Q LN/2 S {p) X iv + i{p,L-x) = Z NpxLf (Q;v) (2.9) 

0<2j+Kp 

where Z^ pX l f (Q;v) is the partition function of the Q-state Potts model on a strip of size 
Np x Lp with cylindrical boundary conditions, as computed in Refs. [39,40]. We have made 
this check for p = 4, 5, 6 and for several values of L and N. 

2) For p = 3 the RSOS model trivializes. Only the xi,i sector exists, and T x is one- 
dimensional for all L. Eq. (2.5) gives simply 

Z LfxNp (Q = 1;x) = (1+x) e , (2.10) 

where E is the number of lattice edges (faces on the height lattice). It is not possible to treat 
the bond percolation problem in the RSOS context, since this necessitates taking Q — > 1 as 
a limit, and not to sit directly at Q — 1. Hence, the right representation for studying bond 
percolation is the FK representation. 

3 Square- lattice Potts model with free cyclic boundary 
conditions 

3.1 Ising model (p = 4) 

The partition function for a strip of size L F x Np is given in the RSOS representation as 
Z LfxNp (Q = 2 ;x) = 2 NL ' 2 \ Xlyl {x) + Xl ^x)] (3.1) 

7 For exceptional values of x there may still be cancellations between eigenvalues with opposite sign. 
However, the pair of eigenvalues that cancel must now necessarily belong to the same sector xi,2j+i- 
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where Xi,2j+i( x ) = trT 2 j+i(p = 4:,L;x) N . The dimensionality of the transfer matrices can 
be obtained from the general formulae derived in Appendix A: 

dimT fc (p = 4,L) = 2 L ~ X , A; = 1,3 (3.2) 

We have computed the limiting curves B L for L = 2,3,4. These curves are displayed 
in Figure 5(a)-(c). In Figure 5(d), we show simultaneously all three curves for comparison. 
In addition, we have computed the partition-function zeros for finite strips of dimensions 
Lp x (pL)p for aspect rations p = 10, 20, 30. These zeros are also displayed in Figure 5(a)- 
(c). 8 For 5 < L < 8, we have only computed selected features of the corresponding limiting 
curves (e.g., the phase diagram for real x). 



3.1.1 L = 2 

This strip is displayed in Figure 3. Let us denote the basis in the height space as 
\h\,h,2,h^hi,h^), where the order is given as in Figure 3. 

The transfer matrix Ti is two-dimensional: in the basis {|1, 2, 1, 2, 1), |1, 2, 3, 2, 1)}, it 
takes the form 

T ><*=^ = 2 > = ^(%£) < 3 - 3 ' 

where we have used the shorthand notation 

Y L<k = x k (x + V2 j , k = 0, . . . , 2L - 1 (3.4) 

The transfer matrix T 3 is also two-dimensional: in the basis {|1, 2, 1, 2, 3), |1, 2, 3, 2, 3)}, 
it takes the form 

™ = 4,* = *> = ^(%£) (3-5) 
For real x, there is a single phase-transition point at 

x c = — l — « -0.7071067812 (3.6) 
v2 

This point is actually a multiple point. 9 There is an additional pair of complex conjugate 
multiple points at x — — e ±Mr / 4 = — l/y/2 ± i/y/2. We also find an isolated limiting point at 
x = —\/2 due to the vanishing of all the eigenvalues (see Ref. [31] for an explanation of this 
issue in terms of the Beraha-Kahane- Weiss theorem). 

The dominant sector on the real x-axis is always except at x — —\pl and x — — 1/ \/2; 
at these points the dominant eigenvalues coming from each sector xi,k become equimodular. 
On the regions with null intersection with the real x-axis, the dominant eigenvalue comes 
from the sector Xi,3- 

8 After the completion of this work, we learned that Chang and Shrock had obtained the limiting curves 
for L = 2 [41, Figure 20] and L = 3 [42, Figure 7]. The eigenvalues and amplitudes for L = 2 had been 
previoulsy published by Shrock [43, Section 6.13]. 

9 Throughout this paper a point on Bl of order > 4 is referred to as a multiple point. 
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3.1.2 L > 3 



For 3 < L < 8, we find two phase-transition points on the real a;- axis: 

1 

x c ,i = -= « -0.7071067812 (3.7a) 

v2 

Xc 2 = -y/2 « -1.4142135624 (3.7b) 

Both points are actually multiple points (except x Cj 2 for L = 3). There is an additional pair 
of complex conjugate multiple points at x — — e ±M "/ 4 . 

For x > x Ct i, the dominant eigenvalue always belongs to the sector xi,i- For x < x C) i, 
this property is true only for even L = 4,6, 8; for odd L = 3, 5, 7, the dominant eigenvalue 
for x < x Ct i belongs to the Xi,3 sector. 



3.2 Q = B 5 model (p = 5) 

The partition function for a strip of size L F x N P is given in the RSOS representation as 



Z Lf xN p (Q = B 5 ;x) = B^ L/2 Xi,i( x ) + V B 5Xi,3( x ) 



(3. 



where Xi,2j+i(x) = tr J 2 j+i(p = 5, L; x) N . 

We have computed the limiting curves Bl for L = 2, 3, 4. These curves are displayed in 
Figure 6(a)-(c). In Figure 6(d), we show all three curves for comparison. For L = 5,6, we 
have only computed selected features of the corresponding limiting curves. 



3.2.1 L = 2 



The transfer matrix T x is two-dimensional: in the basis {|1, 2, 1, 2, 1), |1, 2, 3, 2, 1)}, it 
takes the form 



T!(p = 5,L = 2) 



^/B*X 3 Bl /4 ^B*xXl 
Bl 1//4 x 3 x 2 (l + x) 



where we have used the shorthand notation 

X 3 = x+y^Ek, X* 3 = x+^/B* 
in terms of B 5 and B£ defined as 

3 + V5 3-^5 



Bt 



(3.9) 



(3.10) 



(3.11) 



The transfer matrix T 3 is three-dimensional. In the basis {|1, 2, 1, 2, 3), |1, 2, 3, 4, 3), 
|1, 2, 3, 2, 3)}, it takes the form 



T 3 (p = 5,L = 2) 



/ ,fB% xX\ 
fBjx 2 



o b; 1/4 x 2 x 3 \ 

xX* B* 5 1/4 x(l+x) 
\Bl lfA x 2 (l + x) Bl 1/A x x{l + xf ) 



(3.12) 



For real x, there is a single phase-transition point at 

VB~ 5 1 + V5 
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-0.8090169944 



(3.13) 



We have also found that the limiting curve contains a horizontal line between x = xbk = — 1 
and x ~ —1.3843760945. The latter point is a T point, and the former one, a multiple point. 
There is an additional pair of complex conjugate multiple points at 

x = - e ±iw / 5 = ± 1/4 « -0.8090169944 ± 0.5877852523 i (3.14) 

We have found two additional pairs of complex conjugate T points at x ~ —1.5613823329 ± 
0.3695426938 i, and a; « -0.9270509831 ± 0.3749352940 %. The dominant sectors on the real 
x-axis are 

• xi,i for x G (-oo, -1.3843760945) U (-0.8090169944, oo) 

• xi,3 for a; G (-1.3843760945, -0.8090169944) 



3.2.2 L > 3 

For L = 3 there are two real phase-transition points at 

x cA « -2.1862990086 (3.15a) 
x Ci 2 « -0.9176152641 (3.15b) 

The limiting curve contains a horizontal line between two real T points x ~ —1.2066212246 
and x ks —0.9713270390. There are nine additional pairs of complex conjugate T points. 
The dominant sectors on the real x-axis are 

• xi,i for x e (-2.1862990086, -1.2066212246) U (-0.9176152641, oo) 

• xi,3 for x G (-oo, -2.1862990086) U (-1.2066212246, -0.9176152641) 
For L — 4, the real transition points are located at 

x c> i ss -1.3829734471 (3.16a) 
x c ,2 ~ -0.9475070976 (3.16b) 

We have found that the curve B4 contains a horizontal line between two real T points: 
x ~ —1.1982787848 and x ~ —0.9776507663. Two points belonging to such line are actually 
multiple points: x w -0.9923357481 and a: w -0.9972135728. We have found 34 pairs of 
complex conjugate T points. The phase diagram is rather involved, and we find several tiny 
closed regions. The dominant sectors on the real x-axis are 

• xi,i forxG (-00, x cA ) U (-1.1982787848, -0.9972135728) U (x c>2 , 00) 

• xifl for x G (x c ,i, -1.1982787848) U (-0.9972135728, x c>2 ) 
For L = 5, there are four real phase-transition points at 



x Cyl w -2.4492425881 (3.17a) 

x Ci2 w -1.2097913730 (3.17b) 

x Ci3 w -1.1717714277 (3.17c) 

a; Ci4 w -0.9616402644 (3.17d) 



Again, i3 5 contains a horizontal line between x w —1.1323655119 and x ps —0.9770339631. 
The dominant sectors on the real x-axis are 
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• xi,i for x e (^c,n -0.9770339631) U (ar C)4 , oo) 

• xi,3 for x G (-oo,x C) i) U (-0.9770339631, x c>4 ) 

Finally, for L = 6, there are five real phase-transition points at 



x C) i « -1.2750054535 (3.18a) 

x Ci 2 « -1.2712112920 (3.18b) 

x Ci 3 « -1.1323753929 (3.18c) 

x c>4 « -1.1052066740 (3.18d) 

x c<5 w -0.9700021428 (3.18e) 



The limiting curve contains a horizontal line between two real T points: x pa —1.0877465961 
and x pa -0.9792223546. This line contains the multiple point x pa -1.0781213888. The 
dominant sectors on the real x-axis are 

• xi,i for x e (-oo, x c ,i) U (x c ,2, -1.0877465961) U (-1.0781213888, oo) 

• xi,3 for a: G (x c ,i, x c , 2 ) U (-1.0877465961,-1.0781213888) 

In all cases 2 < L < 6, there is a pair of complex conjugate multiple points at x = 
_ e ±«r/5 ^ -0.8090169944 ± 0.5877852523 i. 



3.3 Three-state Potts model (p = 6) 

The partition function for a strip of size L F x iV P is given in the RSOS representation as 

Z LfxNp (Q = 3;x) = 3 NL /*[ Xl!l (x)+2 X i,s(x)+XiAx)} (3-19) 

where Xi,2j+i(^) = tr T 2 j+i(p = 6, L; x) N . 

We have computed the limiting curves Bl for L = 2,3,4. These curves are displayed in 
Figure 7(a)-(c). 10 

In Figure 7(d), we show all three curves for comparison. For L = 5, 6, 7 we have only 
computed selected features of the corresponding limiting curves. 

3.3.1 L = 2 

The transfer matrix T 5 is one-dimensional, as there is a single basis vector {|1, 2, 3, 4, 5)}. 
The matrix is given by 

T 5 (p = 6,L = 2) = x 2 (3.20) 

The transfer matrix Ti is two-dimensional: in the basis {|1, 2, 1, 2, 1), |1, 2, 3, 2, 1)}, it 
takes the form 

1 / Xf V2xXf 



10 After the completion of this work, we learned that the limiting curves for the smallest widths had been 
already obtained by Chang and Shrock: namely, L = 2 [42, Figure 22], and L — 3 [42, Figure 8]. Please note 
that in the latter case, they used the variable u = l/(v + 1) = l/(x\/Q + 1), instead of our variable x. 
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where we have used the shorthand notation 

X 1 = x + V3, X 2 = 2x + V3 (3.22) 

The transfer matrix T 3 is three-dimensional. In the basis {|1, 2, 1, 2, 3), |1, 2, 3, 4, 3), 
|1, 2, 3, 2, 3)}, it takes the form 



/ 2xX\ 2 v / 2x 2 X 1 
T 3 (p = 6,L = 2) = :ri = | ^/6x 2 t/3xX 2 V^xX 2 | (3.23) 

y/2x 2 X 2 3a; xXj 



2^3 



For real x, there are two phase-transition points 

Xc>l = -V3 = x_ « -1.7320508076 (3.24a) 
a/3 

x C)2 = « -0.8660254038 (3.24b) 

There is one pair of complex conjugate T points at x « —1.6522167507 ± 0.5104474197i. 
There are three multiple points at x — —y/3/2, and x = — V3/2 ± i/2 = — e ±l7r / 6 . The 
dominant sectors on the real x-axis are 

• xi,i f° r x £ ( — °°, —V%) U (—a/3/2, oo) 

• xi,5 for x G (->/3, -V3/2) 

On the regions with null intersection with the real x-axis, the dominant eigenvalue comes 
from the sector Xi 3- 



3.3.2 L > 3 

For L = 3, there are three real phase-transition points 

x C;1 « -1.9904900679 (3.25a) 
a; C)2 = _V3 = x_ « -1.7320508076 (3.25b) 

x C)3 = « -0.8660254038 (3.25c) 

The limiting curve contains a small horizontal segment running from x ~ —1.0539518478 
to x = xbk = — 1- On this line, the two dominant equimodular eigenvalues come from the 
sector xi,5- 

We have found 15 T points (one real point and seven pairs of complex conjugate T 
points). The real point is x = —1. The phase structure is vastly more complicated than that 
for L = 2. In particular, it contains three non-connected pieces, and four bulb-like regions. 
On the real x-axis, the dominant eigenvalue comes from 

• xi,i f° r x e (x Cj i, —\/3) U (-V3/2, oo) 

• xi,3 for x e (-oo,x C) i) U (-^,-1.0539518478) 

• xi,s for x e (-1.0539518478, — s/3/2) 
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For L — 4, there are four phase-transition points 



x c ,i = -v 7 ^ = x_ ps -1.7320508076 (3.26a) 

x c ,2 ~ -1.3678583305 (3.26b) 

x C)3 ps -1.2237725061 (3.26c) 
a/3 

x c ,4 = -- y w -0.8660254038 (3.26d) 



This is the strip with smallest width for which a (complex conjugate) pair of endpoints 
appears: x ~ — 0.9951436066±0.00444309186i. These points are very close to the transition 
point x B k = — 1- We have found 36 pairs of conjugate T points. We have also found three 
multiple points at x — —V3, and x = —y/3/2±i/2. The dominant sectors on the real x-axis 
are 

• xi,i for x e (-oo, z Cj3 ) U (—a/3/2, oo) 

• xi,5 for x G (x c> 3, -y/3/2) 

For L = 5, there are six real phase-transition points 



x c ,i w -2.3018586529 (3.27a) 

a; Ci2 = -V3 = x_ ps -1.7320508076 (3.27b) 

x C)3 ps -1.4373407728 (3.27c) 

x cA ps -1.3412360954 (3.27d) 

x C)5 w -1.2613579653 (3.27e) 
/3 

a; Cj6 = ps -0.8660254038 (3.27f) 



We have also found a horizontal line running between the T points x ~ —1.0226306002 and 
x ps —0.9984031794. The dominant sectors on the real x-axis are 

• xi,i for x G (x cA ,x Ct3 ) U (-■ \/3/2,oo) 

• xi,3 for x G (-oo,a; C) i) U (x C)3 , -1.0226306002) 

• xi,5 for x G (-1.0226306002, — V3/2) 

For L — 6, there are also six phase-transition points on the real axis 



x c ,i = -V3 — x_ ps -1.7320508076 (3.28a) 

x Ci2 ps -1.2852299467 (3.28b) 

x C)3 ps -1.2238569234 (3.28c) 

x C)4 ps -1.1271443188 (3.28d) 

x C)5 ps -1.0085262838 (3.28e) 

AT 

x Cj6 = ps -0.8660254038 (3.28f) 



The dominant sectors on the real x-axis are 
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• xi,i for x e (-oo,x Cj2 ) U (x C)3 ,x C)5 ) U (-a/3/2,00) 

• xi,5 for x e (x Ci2 ,x Ci3 ) U (x c>5 , -a/3/2) 

In all cases 2 < L < 6, we have found three multiple points at x — —a/3, and x = 
-a/3/2 ± i/2 = - e ±i7r / 6 . 

3.4 Four-state Potts model (p = oo) 

It follows from the RSOS constraint and the fact that ho — 1 is fixed, that the maximal 
height participating in a state is h max = max(2L,p — 1). In particular, for any fixed L the 
number of states stays finite when one takes the limit p — > oo. Meanwhile, the Boltzmann 
weight entering in Eq. (2.7) has the well-defined limit (hjh'j) 1 ^ 2 /hj-i, and the amplitudes 
(2.2) tend to Sj (oo) = 2j + 1. We shall refer to this limit as the p = oo (or Q — 4) model. 

We have computed the limiting curves Bl for L = 2,3,4. These curves are displayed in 
Figure 8(a)-(c). In Figure 8(d), we show all three curves for comparison. 

3.4.1 L = 2 

The transfer matrices are 



1 / (x + 2) 3 VZx(x + 2y 
Tl " 2 *v V3x* x*(2 + 3x) ) (3 ' 29a) 



1 



3a;(a; + 2) 2 3\/3x 2 (:r + 2) 



T 3 = - 2^6x 2 2x(3a; + 4) 2a/2 :r(3x + 2) | (3.29b) 

1 V3x 2 (3x + 2) 4a/2 x x(3x + 2) 2 

T 5 = x 2 (3.29c) 

For real x, we find a multiple point at x = —1, where all eigenvalues become equimodular 
with |Aj| = 1. The dominant sector on the real x-axis is always Xi,i- 

3.4.2 L > 3 

For L = 3 there are two real phase-transition points: x — — 1 (which is a multiple point), 
and x c ~ —1.6424647621. We have found ten pairs of complex conjugate T points and a 
pair of complex conjugate endpoints. The dominant sectors on the real x-axis are Xi,3 f° r 
x < —1, and xi,i f° r x > ~ 1- The sector xi,5 is only dominant in two complex conjugate 
regions off the real x-axis, and the sector xij is never dominant. 

For L = 4 we only find a single real phase-transition point at x — — 1. We have also 
found 32 pairs of complex conjugate T points and two pairs of complex conjugate endpoints. 
The dominant sector on the real x-axis is always Xi,i- There is also two complex conjugate 
regions where the dominant eigenvalue comes from the sector xi,5, an d the sectors Xi,7 an d 
Xi,9 are never dominant in the complex x-plane. 

For L = 5 we find four real phase-transition points at 

Xcl = -1.9465787472 (3.30a) 

x C)2 = -1.5202407889 (3.30b) 

x C)3 = -1.3257163278 (3.30c) 

x C)4 = -1 (3.30d) 
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The dominant sectors are Xi,3 f° r x £ ( — oo,x Cj i) U (^c,2, — 1); and m the region x G 
(x c ,i,x C:2 ) U (-l,oo). 

For L = 4 we only find a single real phase-transition point at re = —1. The dominant 
sector on the real x-axis is always Xi,i- 

In all cases 3 < L < 5, the point x = — 1 is a multiple point where all the eigenvalues are 
equimodular with |Aj| = 1. 

4 Common features of the square-lattice limiting curves 
with free cyclic boundary conditions 

From the numerical data discussed in Sections 3.1-3.3, we can make the following con- 
jecture that states that certain points in the complex x-plane belong to the limiting curve 

Conjecture 4.1 For the square-lattice Q-state Potts model with Q = B p and widths L > 2: 

1. The points x = — e ±l7r//p belong to the limiting curve. At these points, all the eigenvalues 
are equimodular with \\\ = l. 11 Thus, they are in general multiple points. 

2. For even p, the point x = —y/Q/2 always belongs to the limiting curve Bl- 12 Further- 
more, if p = 4, 6, then the point x = —\[Q also belongs to Bl- 

The phase structure for the models considered above show certain regularities on the real 
a;-axis (which contains the physical regime of the model). In particular, we conclude 

Conjecture 4.2 For the square-lattice Q-state Potts model with Q — B p and widths L > 2: 

1. The relevant eigenvalue on the physical line v G [—1, 00) comes from the sector xi,i- 

2. For even L, the leading eigenvalue for real x comes always from the sector Xi,i> except 
perhaps in an interval contained in [— \J~Q, —\fQj2\. 

3. For odd L, the leading eigenvalue for real x comes from the sector xi,3 for all x < x < 
— \[Q, and from the sector Xi,i for all x > —\fQj2. 

In the limiting case p = 00 the RSOS construction simplifies. Namely, the quantum 
group U q (SU(2)) reduces to the classical U(SU(2)) (i.e., q — * 1), and its representations no 
longer couple different K 1:2 j + i, cf. Eq. (2.4). Accordingly we have simply K 1:2 j+i = Xi,2j+i- 
When increasing p along the line xbk(Q), the sector Ki^j+i which dominates for irrational 
p will have higher and higher spin j [7]; this is even true throughout the Berker-Kadanoff 
phase. 13 One would therefore expect that the p = 00 RSOS model will have a dominant 
sector Xi,2j+i with j becoming larger and larger as one approaches xbk(Q = 4) = —1. 

11 This property has been explicitly checked for all the widths reported in this paper. 
12 This property has been verified for p = 8, 10 and 2 < L < 6. 

13 See Ref. [34] for numerical evidence along the chromatic line x = —\j\[Q which intersects the BK phase 
up top = 12 [15]. 
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This argument should however be handled with care. Indeed, for p — > oo the BK phase 
contracts to a point, (Q,v) = (4, —2), and this point turns out to be a very singular limit 
of the Potts model. In particular, one has xbk — x± for Q — 4, and very different results 
indeed are obtained depending on whether one approaches (Q, v) = (4, —2) along the AF 
or the BK curves (1.4). This is visible, for instance, on the level of the central charge, with 
c — > 2 in the former and c — > — oo in the latter case. To wit, taking x — > — 1 after having 
fixed p = oo in the RSOS model is yet another limiting prescription, which may lead to 
different results. 

The phase diagrams for Q = 4 (p — > oo) do agree with the above general conjectures 4.1- 
4.2. In particular, whenp — * oo, the multiple points — e ±Mr / p — > — 1 = x B k (Conjecture 4.1.1) 
and this coincides with the point — ^jQ/2 (Conjecture 4.1.2). On the other hand, the sector 
is the dominant one on the physical line v G [—1, oo) (Conjecture 4.2.1), and we observe 
a parity effect on the unphysical regime v G (— oo, —1). For even L, the only dominant sector 
is xi,i m agreement with Conjecture 4.2.2 (although there is no interval inside [—2, —1] where 
Xi,3 becomes dominant). For odd L, Conjecture 4.2.3 also holds with x = —\fQ = —2 (at 
least for L — 3, 5). For L = 2, 3, 4, we find that in addition to the sectors xi,i an d Xi,3> on ly 
the sector xi,5 becomes relevant in some regions in the complex x-plane. 



4.1 Asymptotic behavior for \x\ — > oo 

Figures 5-8 show a rather uncommon scenario: the limiting curves contain outward 
branches. As a matter of fact, these branches extend to infinity (i.e., they are unbounded 14 ), 
in sharp contrast with the bounded limiting curves obtained using free longitudinal boundary 
conditions [39,40]. It is important to remark that this phenomenon also holds in the limit 
p — > oo, as shown in Figure 8. 

As \x\ — > oo these branches converge to rays with definite slopes. More precisely, our 
numerical data suggest the following conjecture: 15 

Conjecture 4.3 For any value of p, the limiting curve Bl for a square-lattice strip has 
exactly 2L outward branches. As \x\ — > oo ; these branches are asymptotically rays with 

argx = 9 n (L) = tt ^ - -L) , n = 1, 2, . . . , 2L (4.1) 

By inspection of Figures 5-8, it is also clear that the only two sectors that are relevant 
in this regime are xi,i an d Xi,3- I n particular, the dominant eigenvalue belongs to the xi,i 
sector for large positive real x, and each time we cross one of these outward branches, the 
dominant eigenvalue switches the sector it comes from. In particular, we conjecture that 

Conjecture 4.4 The dominant eigenvalue for a square-lattice strip of width L in the large 
\x\ regime comes from the sector Xi,i in the asymptotic regions 

argx G (0 2n _!(L),0 2n (L)) , n=l,2,...,L (4.2) 

In the other asymptotic regions the dominant eigenvalue comes from the sector xi,3- 

14 An unbounded branch is one which does not have a finite endpoint. 

15 Chang and Shrock [42] observed for L = 3 that if we plot the limiting curve in the variabe u = 
l/(x^/Q+ 1), then the point u = is approached at specific angles argw consistent with our Conjecture 4.3. 
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In particular, this means that for large positive x the dominant sector is always xi,i- 
However, for large negative x the dominant eigenvalue comes from xi,i is L is even, and 
from Xi,3 if L is odd. Thus, this conjecture is compatible with Conjecture 4.2. 

An empirical explanation of this fact comes from the computation of the asymptotic 
expansion for large \x\ of the leading eigenvalues in each sector. It turns out that there is a 
unique leading eigenvalue in each sector xi,i an d Xi,3 when \x\ — > oo. As there is a unique 
eigenvalue in this regime, we can obtain it by the power method [44]. Our numerical results 
suggest the following conjecture 

Conjecture 4.5 Let K tl (L) (resp. \*^(L)) be the leading eigenvalue of the sector xi,i (resp. 
Xi,s) i n regime \x\ — > oo. Then 



a k (L) k 

Jb 



., t Q k/2 



(4.3a) 



KM = Q^'^x^- 1 

K,i(L) - KM = VQx^ + SiL-l^ + Oix 1 - 3 ) (4.3b) 
Furthermore, we have that 

ai(L) = 2L-1, L>2 (4.4a) 

a 2 (L) = 2L 2 - 3L + 1 , L > 3 (4.4b) 

The first coefficients a fe (L) are displayed in Table 1; the patterns displayed in (4.4) are 
easily verified. The coefficients a k (L) also depend on p for k > 3. 

Indeed, the above conjecture explains easily the observed pattern for the leading sector 
when x is real. But it also explains the observed pattern for all the outward branches. These 
branches are defined by the equimodularity of the two leading eigenvalues 

|A*,i| = |A* >3 | = \K,i-VQx L - l + 0{x L - 2 )\ (4.5) 

This implies that 

Re[A^ix L_1 ] = (4.6) 

where x is the complex conjugate of x. Then, if x = \x\e ld , then the above equation reduces 
to 

cos (9L) = 6 n = ^(2n - 1) , n = 1, 2, . . . , 2L (4.7) 
in agreement with Eq. (4.1). 

Remark. The existence of unbounded outward branches for the limiting curve of the Potts 
model with cyclic boundary conditions is already present for the simplest case L — 1. Here, 
the strip is just the cyclic graph of n vertices C n . Its partition function is given exactly by 

Z Cn (Q,v) = (Q + v) n + (Q-l)v n (4.8) 

Then, we have two eigenvalues Ai = Q + v = Q + x^/Q and K — v — XyfQ, which grow 
like ~ x 2L ~ x = x and whose difference is Q = 0{x L ~ v ), in agreement with Conjecture 4.5. 
Furthermore, the limiting curve is the line Rex = which, as |x| — > oo, has slopes 

given by ±7r/2, in agreement with Conjecture 4.3. 
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4.2 Other asymptotic behaviors 

For the Ising case (p = 4) the points x = —\J~2 and x = — \j\[2 are in general multiple 
points and we observe a pattern similar to the one observed for \x\ — > oo. 

For x = —l/y/2, we find that, if we write x = — l/v^ + e with |e| <C 1, within each sector 
there is only one leading eigenvalue A*j(L) ~ 0(1). More precisely, for L > 3, 

K,i(L) = 2- L / 2 + 0(e 3 ) (4.9a) 
A* fl (L) - A*, 3 (L) = 2e L + 0{e L+1 ) (4.9b) 

Again, the equimodularity condition when |e| — > implies that Re(e L ) = 0, whence arge = 6 n 
with 9 n given by Eq. (4.7). 

The case x = —y/2 is more involved. If we write x = —\/2 + e with |e| <C 1, we find that 
in the sector xi,i there are two eigenvalues of order 0(e), and the rest are of order at least 
0(e 2 ). The same conclusion is obtained from the sector xi,3- If we call A^ (i = 1,2) the 
dominant eigenvalues coming from sector xij, then we find for L > 3 that 

A (i) (L) = 2 ( L -m e + ^) « -A$(L) (4.10a) 

$l(L) + $l(L) = \^\ J LeV6n (4.10b) 
*' 1V ' *' 1V ; [2(L-l)e L L odd V ; 

A«(L) = - 2 W*e + 0(<?) » -A3(L) (4.10c) 

ASW + A«(L) = t^L e L-i + {e L ) (4.10e) 

The equimodularity condition implies that 

Refee 1 " 1 ] = cos(0(L - 2)) = (4.11) 
Thus, the same asymptotic behavior is obtained as for x = —l/y/2, except that L — > L — 2: 

^ n = 2^2) (2n_1) ' n=l,...,2(L-2) (4.12) 

5 Triangular-lattice Potts model with free cyclic bound- 
ary conditions 

5.1 Ising model (p = 4) 

For this model we know [16-18] the exact transition temperature for the antiferromagnetic 
model i> c ,af = — 1 = v c ,bk- The partition function is given by a formula similar to that of 
the square lattice, and the dimensionality of Tj(2,L) is the same as for the square lattice. 
In what follows we give the different matrices in the same bases as for the square lattice. 

We have computed the limiting curves Bl for L = 2,3,4. These curves are displayed in 
Figure 9(a)-(c). 16 In Figure 9(d), we show all three curves for comparison. 

16 After the completion of this work, we learned that Chang and Shrock had obtained the limiting curve 
for L = 2 [41, Figure 18]. 
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5.1.1 L = 2 



This strip is drawn in Figure 4. The transfer matrices are 

1 / 2x 4 + 5V2x 3 + 12x 2 + 8y/2x + 4 x(2x 3 + hV^x 2 + 8x + 2y/2) i , , 
1 ~ 2 V x 2 (2x 2 + 3 v / 2x + 2) x 2 (2x 2 + 3v^x + 2) - 1 



2x 3 + Sv^x 2 + 8x + 2^ 2x 3 + 5V2x 2 + 8x + 2y/2 \ 



Ts 2 V x{2x 2 + ?>V2x + 2) 8x 3 + Sv^x 2 + 6x + 2^ ) ( °' J1)) 
For real x, there is a single phase-transition point at 

x c = -l/y/2 « -0.7071067812 (5.2) 
We have found that the entire line 

Rex = -l/y/2 (5.3) 

belongs to the limiting curve. Furthermore, B2 is symmetric with respect to this line. Finally, 
there are two complex conjugate multiple points at x = — l/y/2 ± i/y/2 = — e ±Mr / 4 . 

The dominant sector on the real x-axis is xi,i for x > — 1/\/2, and Xi,3 for x < — l/y/2. 
Note that x c = -l/S/2 gives the right bulk critical temperature for this model in the anti- 
ferromagnetic regime. 



5.1.2 L > 3 

For L = 3,4 we have found that a) The line Rex = — l/y/2 belongs to the limiting 
curve; b) Bl is symmetric under reflection with respect to that line; c) Bl contains a pair 
of multiple points at x = — e ±Mr / 4 ; and d) The dominant sector on the real x-axis is Xi,i f° r 
x > — l/y/2, and Xi,3 f° r x < — 1/>/2. 

For L = 3, there is another pair of multiple points at x « -l/v 7 ^ ± 0.7257238112 i; for 
L = 4 this pair is located at x w -l/y/2 ± 0.7647261156 i. 

For L = 5, 6, 7, we have found that there is a single real phase-transition point at x = 
— l/y/2, and that the dominant sector for x > — l/y/2 (resp. x < —l/y/2) is Xi,i (resp. Xm)- 



5.2 Q = B 5 model (p = 5) 

We have computed the limiting curves B L for L = 2,3,4. These curves are displayed in 
Figure 10(a)-(c). In Figure 10(d), we show all three curves for comparison. 



5.2.1 L = 2 

The transfer matrices are 



Ti 



x 



B 5 + x(2x + 4X 3 + x 2 Xl) xB, 



x'Bl^il + Sy/Bfx + x 2 ) 



f B~ b + Ax + x 2 Xl 



1/4/ 



5 VVJ B 5 + 4x + x 2 X*) 
x 2 (l + 3x + 3 v / ^x 2 ) 

Bl /4 (l + xy/5B 5 



(5.4a) 



>^x 2 Xl) 

{Bty'\i + x) 

(B*) 1 ^! + 3x + y/B~ 5 x 2 ) (B*y/ 4 (l + x) 1 + 3x + 3x 2 + y/B~ 5 x 3 



'B%x 



(5.4b) 
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where we have defined the shorthand notations 

XI = l + 3y/B* + X* (5.5a) 
X* 5 = l+A^Bl + Xt (5.5b) 

For real x, there are two phase-transition points at 

x c>1 « -0.9630466372 (5.6a) 
x c , 2 ~ -0.5908569607 (5.6b) 

In fact both points are T points and the whole interval [x Ci i,x Ci 2] belongs to the limiting 
curve £>2. Finally, there are two complex conjugate multiple points at x — — e ±Mr / 5 , as for 
the square-lattice case. The dominant sector on the real x-axis is xi,i f° r x > x c,i, and xi,3 
for x < x C) i. 

5.2.2 L > 3 

For L = 3, there are two real phase-transition points at 

x Cjl « -1.0976251052 (5.7a) 
x c , 2 ~ -0.6376476917 (5.7b) 

We have found two pairs of complex conjugate endpoints at x « — 0.4297467004±0. 6445268125 i, 
and a; ~ —0.3955590901 ± 0.8536454650 %. There are nine pairs of complex conjugate T 
points, and two complex conjugate multiple points at x — — e ±Mr / 5 . The dominant sectors 
on the real x-axis are for x > x c>1 , and Xi,3 f° r x < x c,i 
For L — 4, there are three real phase-transition points at 

x Cjl w -1.0953543257 (5.8a) 
x c , 2 ~ -0.9708876996 (5.8b) 
x c ,3 w -0.6102005246 (5.8c) 

The points x c ^ and a; Ci 3 are T points, and they define a line belonging to the limiting curve. 
This line contains two multiple points at x w -0.6319374252, and x w -0.7685805289. We 
have found two additional pairs of complex conjugate endpoints at x —0.9270404586 ± 
0.3749352143?, and x = — e ±l7r / 5 . In addition, there are 22 pairs of complex conjugate T 
points. The dominant sectors on the real x-axis are 

• for x £ (^0,1^0,2) U (x Cj3 , 00) 

• xi,3 for a; G (-oo,x Cjl ) U (x Ci2 ,x C)3 ) 

For L = 5, we have found five real phase-transition points at 

x cA w -1.0945337809 (5.9a) 

x c ,2 ~ -1.0615208835 (5.9b) 

x Cj3 w -0.8629689747 (5.9c) 

x c ,4 ~ -0.6393693994 (5.9d) 

x Cj5 w -0.6362471039 (5.9e) 

The dominant sectors on the real x-axis are 
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• xi,i for x £ (Zc,l> X c,2) U (x Cj3 , oo) 

• xi,3 for x G (-oo,x C) i) U (x C)2 ,x C)3 ) 

For L = 6 the amount of memory needed for the computation of the phase diagram on 
the real x-axis is very large, so we have focused on trying to obtain the largest real phase- 
transition point. The result is x C:1 w —0.6221939194 < —X/y/B^. The sector xi,i dominates 
for all x > x C; i] and for x < x c>1 , the sector xi,3 dominates. 

5.3 Three-state Potts model (p = 6) 

For this model we also know that there is a first-order phase transition in the antiferro- 
magnetic regime at [40,45] 

a: c ,AF(g = 3) = -0.563512(14) (5.10) 

We have computed the limiting curves Bl for L = 2,3,4. These curves are displayed in 
Figure ll(a)-(c). 17 In Figure 11(d), we show all three curves for comparison. 

5.3.1 L = 2 

The transfer matrices are 

1 / x 4 + 2y/3x 3 + 6x 2 + AVSx + 3 xV2(x 3 + 2 v / 3x 2 + Ax + y/3) \ , 
1 ~ 2 V x 2 V2(x 2 + V3x + l) x 2 (2x 2 + 2V3x + l) J { } 

2(x 3 + 2V3x 2 + Ax + V3) V2X 1 V2(2x 3 + A^3x 2 + 7x + ^3) 
V2x X 2 X 2 

x\/2(2x 2 + 2V3x + 1) X 2 Ax 3 + Ay/Sx 2 + 6x + V3 

(5.11b) 

T 5 = x 2 (5.11c) 

For real x, there are two phase-transition points at 

x C)1 = -2/V3 « -1.1547005384 (5.12a) 
x c>2 = « -0.5773502692 (5.12b) 

The latter one is actually a multiple point. There are also a pair of complex conjugate 
multiple points at x — —e ±tn / 6 = — a/3/2 ±«/2. The dominant sectors on the real x-axis are: 
Xi,i for x > -l/y/3, Xi,3 for x < -2/y/Z, and Xi,5 for x E {-2/y/Z, -l/y/3). 

5.3.2 L > 3 

For L = 3, there are three real phase-transition points at 

Xcl = -2/y/S « -1.1547005384 (5.13a) 
x c , 2 ~ -0.9712924104 (5.13b) 
Xc3 = -1/V3 w -0.5773502692 (5.13c) 

17 After the completion of this work, we learned that Chang and Shrock had obtained the limiting curve 
for L = 2 [41, Figure 19]. 
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The latter one is actually a multiple point. We have found two pairs of complex conjugate 
endpoints at x « -0.3495004588±0.6911735024i, and x « -0.2862942369±0.8514701201 i. 
There are 16 pairs of complex conjugate T points. The dominant sectors on the real x-axis 
are Xi,i for x > -l/y/S, Xi,3 for x < -2/y/Z, and xi,5 for x G (—2/ >/3, — 1/>/3). 
For L = 4, there are five real phase-transition points at 

x c>1 = -2/V3 « -1.1547005384 (5.14a) 

x c ,2 ~ -1.0219801955 (5.14b) 

a; Ci3 ss -1.0041094453 (5.14c) 

x cA « -0.7664034488 (5.14d) 

a; Ci5 = -1/V3 « -0.5773502692 (5.14e) 

The points x Ci3 and rr c 4 are T points, while multiple point. We have found a pair 

of complex conjugate endpoints at x ~ —0.3857232364 ± 0.6652216322 i. In addition, there 
are 14 pairs of complex conjugate T points. The dominant sectors on the real a;- axis are 

• xi,i for x > x cA 

• xi,3 for x < -2/ and x G (x C)2 , x C)3 ) 

• xi,5 for x e (-2/y/3,x c , 2 ) U (x Ci3 ,x C)4 ) 

For L — 5, there are five real phase-transition points at 

Xcl = « -1.1547005384 (5.15a) 

x c ,2 « -0.9326923327 (5.15b) 

x c ,3 ~ -0.7350208125 (5.15c) 

x Ci 4 « -0.6186679617 (5.15d) 

Xc5 = -1/V3 « -0.5773502692 (5.15e) 

The dominant sectors on the real x-axis are 

• xi,i for x e (x C)2 , x C)3 ) U (x Cj4 , oo) 

• xi,3 for x < x c ,2 

• xi,5 for x E {x c>3 ,x cA ) 

For L — 6, there are three real phase-transition points at 

x c ,i = « -1.1547005384 (5.16a) 

x c , 2 « -1.0504774228 (5.16b) 
a; Ci 3 = -1/V3 « -0.5773502692 (5.16c) 

We have also found a small horizontal line belonging to the limiting curve Bq and bounded 
by the T points 

x cA « -0.7688389273 (5.17a) 
x c>5 « -0.7646464215 (5.17b) 

The dominant sectors on the real x-axis are 
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• xi,i for x e (-I/a/3, oo) U (x C)4 , x Cj5 ) 

• xi,3 for x G (-00, -2/ v 7 ^) U (x C;2 , x cA ) 

• Xi,5 for a: G (-2/ >/3, x c , 2 ) U (x Cj5 , -1/ >/3) 

In all cases 3 < L < 6, there is a pair of complex conjugate multiple points at x — — e ±l7r / 6 . 

5.4 Four-state Potts model (p = 00) 

We have computed the limiting curves Bl for L = 2,3,4. These curves are displayed in 
Figure 12(a)-(c). In Figure 12(d), we show all three curves for comparison. 

5.4.1 L = 2 

The transfer matrices are 

T 1 ^X 8 (2x 3 + 3x 2 + 6x + 4) V3xX 8 X 7 \ 

1 = 2 I VZx*X 7 x 2 X 6 J (5 ' 18a) 



1 



3xX 8 (2x 2 + 3x + 2) 2-n/6xX 8 >/3a;X 8 X, 
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T 3 = - 2\/6a; 2 2x(4 + 3x) 2^0^(2 + 3x) | (5.18b) 

° VZx 2 X 6 2V2x(3x + 2) xX 9 

T 5 = x 2 (5.18c) 

where we have defined the short-hand notations 

X 6 = 6x 2 + 9x + 2 (5.19a) 

X 7 = 2x 2 + 3x + 2 (5.19b) 

X 8 = x + 2 (5.19c) 

X 9 = lSx 3 + 27x 2 + 18x + 4 (5.19d) 

For real x, we find a multiple point at x = —1, and a T point at x c as —0.5808613334. 
The limiting curve B2 contains the real interval [— l,x c ]. At x = —1, all eigenvalues become 
equimodular with |A«| = 1. 

We have found two additional pairs of complex conjugate T points atx« — 0.9882427690± 
0.0896233991i, and x w -3/4 ± 0.6614378278 i. The dominant sectors on the real x-axis are 
Xi,i for x > x c , and Xi,3 f° r x < x c- We have found no region in the complex x-plane where 
the sector xi 5 is dominant. 



5.4.2 L > 3 

For L = 3 there are two real phase-transition points: x = — 1 (which is a multiple point), 
and x c w —0.8953488450. The limiting curve contains two connected pieces, two pairs of 
complex conjugate endpoints, 12 complex conjugate T points, and one additional pair of 
complex conjugate multiple points at x m —3/4 ± 0.6614378278 %. The dominant sectors on 
the real x-axis are Xi,3 for x < —1; Xi,5 for x G (—1, x c ); and Xi,i for x > x c . We have found 
no region where the sector Xi,7 is dominant. 

For L = 4 there are two real phase-transition points at x = —1 and x = x c ~ 
—0.7107999762, which is a T point. The real line [— l,x c ] belongs to the limiting curve. 
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The dominant sectors on the real x-axis are: Xi,3 f° r x < — lj Xi,7 f° r x e ( — Ij^c); an d 
for a; > x c . We have found a few small regions with dominant eigenvalue coming from the 
sector Xi,5\ but we have found no region where the sector xi,9 is dominant. 

For L = 5 there are again two real phase-transition points at x — — 1 and x = x c ~ 
—0.8004698444, which is a T point. The real line [— l,x c ] belongs to the limiting curve. The 
dominant sectors on the real x-axis are: Xi,3 f° r x < — lj Xi,9 f° r x £ ( — Ij^c); an d f° r 

For L = 6 there are two real phase-transition points at i = —1 and x = x c ~ 
—0.7033434642, which is a T point. The real line [— l,rr c ] belongs to the limiting curve. 
The dominant sectors on the real x-axis are: Xi,3 f° r x < — lj Xi,u f° r x e ( — 1, a^ c ) j an d Xi,i 
for x > x c . 

In all cases, the point x — —1 is a multiple point where all the eigenvalues are equimodular 
with |Aj| = 1. 

6 Common features of the triangular-lattice limiting 
curves with free cyclic boundary conditions 

The results discussed in Sections 5.1-5.3 allow us to make the following conjecture (in 
the same spirit as Conjecture 4.1 for the square-lattice case) that states that certain points 
in the complex x-plane belong to the limiting curve Bl- 

Conjecture 6.1 For the triangular-lattice Q -state Potts model with Q = B p and width 
L > 2: 

1. The points x = —e ±tn / p belong to the limiting curve. At these points, all the eigenvalues 
are equimodular with |Aj| = 1. Thus, they are in general multiple points. 

2. For even p > 6, the point x = always belongs to the limiting curve Bl- 18 
Furthermore, if p = 4,6, then the point x = —1/t/Q also belongs to Bl- 

The phase diagram on the real x-axis (which contains the physical regime of the model) 
shows certain regularities that allow us to make the following conjecture: 

Conjecture 6.2 For the triangular-lattice Q -state Potts model with Q = B p and width 
L > 2: 

1. For even p, the relevant eigenvalue on the physical line v e [— l,oo) comes from the 
sector xi,i- For odd p, the same conclusion holds for all L > L . 19 

2. The relevant eigenvalue belongs to the sector xi,3 for all real x < —2/^/Q. 

The above conjectures also apply to the limiting case p — > oo (i.e., Q — 4). As for 
the square-lattice case, the multiple points — e ±l7T ^ p — > — 1 as p — > oo (Conjecture 6.1.1) in 
agreement with the fact that x = — 1 is a multiple point for Q = 4. Furthermore, this is also 

18 This property has been verified for p = 6 and 2 < L < 7, and for p = 8, 10 and 2 < L < 5. 

19 For p = 5, we find that Lq = 5. For L = 2,4, the relevant eigenvalue belongs to the sector \i,3 on a 
small portion of the antiferromagnetic physical line we [— l,uo]- 
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in agreement with Conjecture 6.1.2, as in this limit, —2/y/Q = —1. The dominant sectors 
for p — > oo also agree with Conjecture 6.2: on the physical line v G [—1, oo) the dominant 
sector is and for x < —1, the dominant sector is Xi,3- More precisely, we can state the 
following conjecture based on the empirical observations reported above: 

Conjecture 6.3 For the triangular-lattice A-state Potts model defined on a semi-infinite 
strip of width L > 2, there exists some x c (L) > —1 such that xi,i is dominant for x > x c (L), 
Xi,2L-i is dominant for —l<x< x c {L), xi,3 is dominant for x < —1. 



6.1 Asymptotic behavior for \x\ — ► oo 

Figures 9-12 show a similar scenario to the one discussed in Section 4: There are several 
unbounded outward branches with a clear asymptotic behavior for large \x\. Again, this 
scenario also holds in the limit p — > oo (See Figure 12). However there are quantitative 
differences with the scenario found for the square lattice. We should modify Conjecture 4.5 
as follows: 

Conjecture 6.4 Let A^i(L) (resp. \^(L)) be the leading eigenvalue of the sector Xi,i (resp. 
Xi,3) in the regime \x\ — > oo. Then 



i+jr b Mx- k 



Q k/2 



(6.1a) 



K,i{L) = Q L -'x 3L - 2 

X^(L)-X^(L) = 2 L - 1 ^Qx L - 1 + (L-l)2 L - 1 x L - 2 + 0(x L ^) (6.1b) 
Furthermore, we have that 

&i(L) = 3L-2, L>2 (6.2a) 

b 2 {L) = ^L 2 -^-L + 3, L>2 (6.2b) 

b 3 (L) = -L 3 — —L 2 + 13L — 4 , L>3 (6.2c) 

The first coefficients b k (L) are displayed in Table 2; the patterns displayed in (6.2) are 
easily verified. The coefficients bk{L) also depend on p for k > 4. 

Conjecture 6.4 explains the number of outward branches in the triangular-lattice case, 
as well as the observed pattern for the outward branches. Again, these branches are defined 
by the equimodularity of the two leading eigenvalues 

|A*,i| = |A*, 3 | = |A*,i - const, x 1 ' 1 + 0(x L ~ 2 )\ (6.3) 

This implies that 

Re [A*^- 1 ] = (6.4) 
Then, if x = \x\e td , the above equation reduces to 

cos(0(2L-l)) = e n = ^ (2/1-1), n=l,2,...,2(2L-l) (6.5) 

Thus, we get the same asymptotic behavior as for the square lattice with the replacement 
L->2L-1. 
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7 Discussion of the results with free cyclic boundary 
conditions 



The results obtained give indications on the phase diagram of the Potts model, as the 
accumulating points of the zeros of the partition function correspond to singularities of the 
free energy. 

Extrapolating the curves obtained to L — > oo in not an easy matter, given that we have 
only access to relatively small L. However, in Sections 4 and 6 we have noted a number of 
features which hold for all L considered, and hence presumably for all finite L and also in 
the thermodynamic limit. 



7.1 Ising model 

The most transparent case is that of the Ising model (p = 4) on the square lattice. Let 
D(x, r) denote the disk centered in x and of radius r. There are then four different domains 
of interest: 

Pi = D(0,1)\D(- V / 2,1) (7.1a) 
V 2 = D(0, 1) n D(- A 1) (7.1b) 
V 3 = D(-V2,1)\D(0,1) (7.1c) 

= C\ (£(0,1) UD(-v / 2,l)) (7.1d) 



V. 
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The L x N strips with even N are bipartite, whence the Ising model possesses the exact 
gauge symmetry J — > —J (change the sign of the spins on the even sublattice). Since the 
limit N — > oo can be taken through even N only, the limiting curves Bl should be gauge 
invariant. In terms of x the gauge transformation reads 

x^ (7.2) 

l + xy/2 

Note that it exchanges T>2 <-> T>4, while leaving T>\ and T> 3 invariant. In particular, the 
structures of Bl around x = —l/y/2 and \x\ = oo discussed in Section 4 are equivalent. 

On the other hand, the duality transformation x — > 1/x is not a symmetry of Bl' this is 
due to the fact that the boundary conditions prevent the lattice from being selfdual. Note 
that the duality exchanges £>i <-> V A and V 2 <-> T> 3 . But whilst there are many branches of 
Bl in T> A) there are none in T>\. 

The Ising model being very simple, we do however expect the fixed point structure on 
the real a;-axis to satisfy duality. Combining the gauge and duality transformations one can 
connect all critical fixed points: 

gauge duality gauge 

^FM — > x+ — > X- — > XBK, (7.3) 

and the first and the last points in the series are selfdual. In the same way, all the non-critical 
(trivial) fixed points are connected: 

x = o | X | = oo ^ x = -.i/yft d ^l y x = -V2, (7.4) 

and the first and the last points in the series are gauge invariant. 
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The reason that we discuss these well-known facts in detail is that the square-lattice 
Ising model is really the simplest example of how taking p rational (here, in fact, integer) 
profoundly modifies and enriches the fixed/critical point structure of the Potts model, as 
compared to the generic case of p irrational. Taking the limit p — > 4 through irrational values 
we would have had three equivalent c = 1/2 critical points, RG repulsive in x, situated at 
xfm and x±; one c = —25/2 critical point, RG attractive in x, situated at xbk] and two non- 
critical (trivial) fixed points, RG attractive in x, situated at x = and \x\ = oo. This makes 
up for a phase diagram on the real x-axis which is consistent in terms of renormalization 
group flows (see the top part of Fig. 2). 

Conversely, sitting directly at p = 4 replaces this structure by the four repulsive c = 1/2 
critical points (7.3) and the four attractive non-critical fixed points (7.4). This again gives a 
consistent scenario, in which notably the BK phase has disappeared (see the bottom part of 
Fig. 2). In other cases than the Ising model (p > 4 integer) we could expect the emergence of 
even more new (as compared to the case of irrational p) fixed points (critical or non-critical), 
which will in general be inequivalent (due in particular to the absence of the Ising gauge 
symmetry). 

Going back to the case of complex x we can now conjecture: 

Conjecture 7.1 Let T> 1 be the domain defined in Eq. (7. Id). Then 

• The points x such that 

Zl f xn p (Q = 2;x) = (square lattice) (7.5) 
for some L and N are dense in C \ T>\. 

• There are no such points in T>\ . 

We now turn to the Ising model on the triangular lattice. We first note that all the limiting 
curves Bl are symmetric under the combined transformation x <-> —x — a/2 and xi,i ^ Xi,3- 
On the level of the coupling constant this can also be written exp(J) — > — exp(J). 

We also conjecture that 

Conjecture 7.2 Let V tri be the interior of the ellipse 

(Re x + 1/V2) 2 + 3 (Im xf = 3/2 . (7.6) 

Then 

• The points x such that 

Zl f xn p (Q = 2;x) = (triangular lattice) (7.7) 
for some L and N are dense in C \ T> tri . 

• There are no such points in £> t ri- 
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7.2 Models with p > 4 

For square-lattice models with p > 4 the phase diagram in the thermodynamic limit 
is expected to be more complicated. We can nevertheless conjecture that the four values 
x c given by Eq. (1.4), and denoted by solid squares in the figures, correspond to phase 
transition points even for Q = B p a Beraha number. Accordingly, these points are expected 
to be accumulation points for the limiting curves Bl, when L — > oo. 

But these four values of x are not the only fixed points. There is a complex fixed point 
structure between X-(Q) and xbk(Q), and between xbk(Q) and x + (Q). This is because for 
Q equal to a Beraha number, the thermal operator is repulsive at xbk{Q) (and not attractive 
as it would have been in the BK phase for irrational p), whereas it remains repulsive at x_ (Q) 
and x + (Q). Therefore, there must at the very least be one attractive fixed point in each 
of the two intervals mentioned, in order for a consistent phase diagram to emerge. Indeed, 
for p even, there are two new fixed points, one of them being conjectured as —^/Q/2 for all 
even p, and the other being equal to —\[Q only for p = 4 and p — 6. But our results for 
finite L are in favor of an even more complicated structure, involving more new fixed points. 
The structure of the phase diagram for p odd is further complicated by the emergence of 
segments of the real x-axis belonging to Bl- It is however uncertain, whether these segments 
will stay of finite length in the L — > oo limit. 

In the models with p = 5, 6, oo and on both the square and triangular lattices, we have 
found strong numerical evidence to conjecture that the partition-function zeros are dense in 
the whole complex x-plane with the exception of the interior of some domain. The shape of 
this domain depends on both p and the lattice structure; and unlike in the Ising case (p — 4), 
we do not have enough evidence to conjecture its algebraic expression [c.f., Conjectures 7.1 
and 7.2]. For the square lattice and fixed p, the limiting curves Bl seem to approach (from 
the outside) the circles (1.5), especially in the ferromagnetic regime Rex > 0. For the 
triangular lattice and p = oo, the limiting curves in Figure 12 seem to approach the circle 



which goes through the bulk critical points x = — 1 and x = 1/2. 

7.3 The region \x\ > 1 

The emergence of unbounded branches of Bl in the region of \x\ ^> 1 is at first sight 
rather puzzling. Because when \x\ is large enough, we should expect the system to be non- 
critical, and thus be described by a unique leading eigenvalue of the transfer matrix. This is 
at least what happens for the 5- state Potts model on a strip with cylindrical or free boundary 
conditions using the Fortuin-Kasteleyn representation [39,40]. 

One of the main reasons for studying the limiting curves in the first place is that we 
wish to use them to detect the critical points of the models at hand. At a conformally 
invariant critical point there should be an infinite spectrum of transfer-matrix eigenvalues 
I A. 1 > |Ai| > . . . that become degenerate according to [46] |Aj/A | ~ exp(— 2iiXi/L) when 
L — > 00, where Xi are critical exponents. The limiting curves just tell us that the two 
dominant eigenvalues become degenerate, and not even with what finite-size corrections. 
Therefore the fact that a point x (even on the real axis) is an accumulation point of Bl is 
not sufficient for x to be a critical point in the sense of the above scaling behavior. 
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The observed behavior for \x\ 3> 1 just shows that the leading eigenvalues in sectors with 
different boundary conditions {xi,i and Xi.,3) come close. This is most transparent in the Ising 
case, where there is a bijection between RSOS heights and dual spins. It is easily seen that 
Xi,i (resp. Xi ;3 ) corresponds to fixed boundary conditions in the spin representation, with 
all the dual spins on the upper/lower rim being fixed as +/+ (resp. +/—)■ On the other 
hand, within a given sector there should be a finite gap between the leading and next-leading 
eigenvalues, in the region \x\ ^> 1, signaling non-critical behavior. 



7.4 Fixed cyclic boundary conditions 

To avoid the (from the point of view of detecting critical behavior) spurious coexistence 
between two different boundary conditions, we should rather pick boundary conditions that 
break the Zq symmetry of the Q-state Potts model explicitly. We now illustrate this pos- 
sibility by making a particular choice of fixed boundary conditions, which has the double 
advantage of generalizing those for the Ising case (as discussed above) and enabling the cor- 
responding Potts model partition function Zl x xn p (Q] v) to be written as a sum of RSOS 
model partition functions. 

Consider first the Potts model partition function Z on the dual lattice, with spins S + 
and S- on the upper and lower exterior dual sites, and at the dual coupling J. Recall that 
the duality relation reads simply vv = Q. If we impose free boundary conditions on S±, we 
have by the fundamental duality relation [1] 

Q y - E/2 - l x E Z(Q; Q/v) = Z(Q; v), (7.9) 

where E (resp. V) is the total number of lattice edges (resp. direct sites). Note that V = LN, 
and that E = 2V — iV (resp. E = 3V — 2N) for the square (resp. triangular) lattice. We now 
claim that this object with fixed and equal values for S± can again be expressed in terms of 
Ki^j+i, for a generic p. The precise relation reads 

L 

Z LxxNp (Q;v) = Q v ~ E/2 x E Z(Q;Q/v) _ = Q LN/2 f3j(p) K 1>2j+ i(p, L; x), (7.10) 

S+ ~ S - 3=0 

which should be compared with Eq. (2.3). We henceforth refer to Zl x xn p {Q] v) as the 
partition function of the Potts model with fixed cyclic boundary conditions (even though it 
would be more precise to say that it is actually the two exterior dual spins that get fixed). 
The amplitudes read 



ftW = if + (-D , (i-qJ ■ P.") 

Note that for arbitrary values of Q, the partition function Z(Q; Q/v) can be defined 

S+=S- 

by its FK cluster expansion on the dual lattice, by giving a weight Q to clusters that do not 
contain any of the two exterior sites, and a weight 1 to clusters containing at least one of 
two exterior sites. Eq. (7.10) is a special case of a more general relation which will be proved 
and discussed elsewhere. 

Now, for p integer, we would like to express Z LxxNp (Q; v) in terms of the Xi,2j+i(p, L; x) 
as we did in the case of free cyclic boundary conditions. But because of the (— l) j in the 
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expression of (3j(p), we have f3 np+ j{p) = f3j(p) and (3( n +i) P -i-j = cf. Eq. (2.4) for the 
case of Sj{p), only if p is even. For p even, we can express 

L(P"2)/2J 

Xi,2j+i(p, ^;^) (p even) (7.12) 

j"=o 

which should be compared with Eq. (2.5). For p odd, there does not appear to exist an 
expansion of Z LxxNp in terms of Xi,2j+i- 

Note in particular that f3\{p) = for any p. This has the consequence of eliminating 
the xi,3 sector from the partition function, and, as we now shall see, modify the \x\ ^> 1 
behavior of the phase diagram. 

8 Square- lattice Potts model with fixed cyclic bound- 
ary conditions 

The limiting curves £>l with fixed cyclic boundary conditions (see Figs. 13-16) are very 
similar to those obtained in Ref. [39] for the Potts model with fully free boundary conditions. 
On the other hand, we have already seen that the Bl with free cyclic boundary conditions 
are very different. 

Before presenting the results for fixed cyclic boundary conditions in detail we wish to 
explain this similarity. We proceed in two stages. First we present an argument why the 
limiting curves corresponding to just the sector xi,i almost coincide with those for fully free 
boundary conditions. Second, we take into account the effect of adding other sectors Xi,2j+i- 

Let Tfk be the transfer matrix in the FK representation with zero bridges (cf. footnote 
6), and let A« be its eigenvalues. 20 Then one has, with cyclic boundary conditions 

^i,i = trT^ K = EAf. (8.1) 

i 

Due to the coupling of K^ 2 j+i, given by Eq. (2.6), the eigenvalues of T 1 (i.e., the transfer 
matrix that generates cf. Eq. (2.8)) form only a subset of the eigenvalues of T FK . More 
precisely, 

Xi,i = J>Af, ( 8 - 2 ) 

i 

where at = or 1 are independent of x. Note that when L < p — 1, Eq. (2.6) gives simply 
Xi,\ = if 1,1, and so in that case all a, = 1. 

Meanwhile, the partition function of the Potts model with fully free boundary conditions 
is given by [32] 

Z {TCC =(f\T§ K \ i ) = J2^^i (8-3) 

i>l 

where the amplitudes ctj are due to the free longitudinal boundary conditions. Note that 
some of the ce, could vanish identically, and indeed many of them do vanish. For example, in 
the case of the square lattice, the vectors \i) and (/| are symmetric under a reflection with 

20 We label the Aj by letting Ao be the eigenvalue which dominates for x real and positive, and using 
lexicographic ordering [32] for the remaining eigenvalues. 
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respect to the axis of the strip, whence only the Aj corresponding to eigenvectors which are 
symmetric under this reflection will contribute to Z^ ee . 

For x > real and positive, it follows from simple probabilistic arguments that the 
dominant eigenvalue A will reside in the zero-bridge sector K X1 and is not canceled by 
eigenvalues coming from other sectors. Therefore a = 1. On the other hand, the Perron- 
Frobenius theorem and the structure of the vectors \i) and (/| implies that a > 0. We 
conclude that the dominant term in the expansions of Xi,i an d Z^ cc are proportional. By 
analytic continuation the same conclusion holds true in some domain in the complex rr-plane 
containing the positive real half-axis. Moving away from that half-axis, a first level crossing 
will take place when A crosses another eigenvalue Aj. If none of the functions ctj and <Sj 
are identically zero, the corresponding branch of the limiting curve Bl coincides in the two 
cases. Further away from the positive real half-axis other level crossings may take place, and 
the limiting curves remain identical until a level crossing between Xj and A& takes place in 
which either otj = and aj ^ 0, or conversely otj ^ and c\j — 0. When L < p — 1 the only 
possibility is the former one, since all a.i — 1. 

If we now compare the limiting curves of Zf rcc and Zrsos, the latter being defined as 
some linear combination of Xi,2j+i (containing the above argument will be invalidated 
if the first level crossing in Zrsos when moving away from the positive half-axis involves an 
eigenvalue from Xi,2j+i with j > 0. 

With free cyclic boundary conditions, Zrsos contains xi,3- The first level crossing involves 
eigenvalues from xi,i an d Xi,3 (cf- the observed unbounded branches) and is situated very 
"close" [cf. Eqs. (4.7) and (6.5) with n — 1] to the positive real half-axis. Accordingly, the 
limiting curves Bl do not at all resemble those with fully free boundary conditions. On the 
other hand, when Xi,3 is excluded (i.e., in the case of fixed cyclic boundary conditions) the 
first level crossing is between two different eigenvalues from the Xi,i sector (see Figs. 13-16). 

8.1 Ising model (p = 4) 

We have studied the limiting curves given by the sector xi,i in the square-lattice Ising 
case. The results are displayed in Figure 13. It is clear that there are no outward branches, 
as there is a unique dominant eigenvalue in the region \x\ ^> 1. Indeed, this agrees with 
the expected non-critical phase. These curves are very similar to those obtained using the 
Fortuin-Kasteley representation for a square-lattice strip with free boundary conditions [39]. 
In particular, for even L = 2,4 we find that these curves do in fact coincide. However, 
for L = 3 we find disagreements; but only in the region Ret> < — 1. Namely, the complex 
conjugate closed regions defined by the multiple points x = — e -Mr//4 and x = —a/2 (see 
Figure 13b) are replaced by two complex conjugate arcs emerging from x = — e~ l7r / 4 . These 
arcs bifurcate at two complex conjugate T points. 

For L = 2 we find two pairs of complex conjugate endpoints at x ~ —0.5558929703 ± 
0. 1923469388 i, and x ss 0.5558929703 ± 1.60656050121 There is a double endpoint at 

X = -y/2. 

For L = 3 we also find two pairs of complex conjugate endpoints at x ~ — 0.5054436896 ± 
0.1404486742 i, and x « 0.9624601506 ± 1. 16277331801 There is a multiple point at q = 
— a/2, and a pair of complex conjugate multiple points at q — —e~ ±m ^. These multiple 
points also appear in L = 4. 

For L = 4 we find two connected components in the limiting curve. There are two 
pairs of complex conjugate T points at q m —1.1111427356 ± 0.8231882219i, and q m 
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-0.9473515724 ± 0.48947792961 We also find four complex conjugate pairs of endpoints at 
q w -0.6052879436±0.3554255102i, q w -0.4820292937±0. 1111133833 i, q « -0.3346743307± 
1.3000737077 i, and g « 1.0790506924 ± 0.8817674400 1 

8.2 Three-state Potts model (p = 6) 

We have studied the limiting curves given by the sectors xi,i an d Xi,5> °f- Eq. (7.12). 
The results are displayed in Figure 14. We have compared these curves with those obtained 
for a square-lattice strip with free boundary conditions [39]. We find that they agree almost 
perfectly in the region Rex > —1. The only exceptions are the tiny complex conjugate 
branches emerging from the multiple points — e~ l7r / 6 for L = 3,4 and pointing to xbk- The 
differences are in both cases rather small and they are away from the real x-axis. In the 
region Rex < —1, however, the differences between the two boundary conditions are sizeable. 
For free boundary conditions the closed regions tend to disappear, or, at least, to diminish 
in number and size. 

9 Triangular-lattice Potts model with fixed cyclic bound- 
ary conditions 

9.1 Ising model (p = 4) 

We have studied the limiting curves given by the sector xi,i i n the triangular-lattice Ising 
case. The results are displayed in Figure 15, and they are the same than those obtained 
with the Fortuin-Kasteleyn representation [40], with free boundary conditions, for all L. 
Therefore, we see a non-trivial effect of the lattice: for the triangular lattice, the dominant 
eigenvalues always comes from K\\, contrary to the case of the square lattice. 

For L = 2 we find two real endpoints at q — — \/2 and q = — l/y/2, and an additional 
pair of complex conjugate endpoints at x ~ 0.3535533906 ± 0.9354143467 1 At x — — 1 there 
is a crossing between the two branches of the limiting curve. 

For L = 3 we find two real endpoints at q — — y/2 and q = —l/y/2, and four pairs of 
complex conjugate endpoints at x w -1.4346151869 ± 0.9530458628 i, x w 0.5477064083 ± 
0.6206108204 i, x w -l/y/2± 0.4918781633 i, and x w -l/ v / 2±0.9374415716 %. The limiting 
curve contains two complex conjugate vertical lines determined by the latter two pairs of 
endpoints, and a horizontal line determined by the two real endpoints. We have found 
three pairs of complex conjugate T points at x — ly/2 ± 0.5353475100 i, x w — ly/2 ± 
0.7246267519 i, and x m -ly/2 ± 0.8539546894 %. Finally, there is a multiple point at x « 
-0.9681295813. 

For L = 4, we again find a horizontal real line bounded by two real endpoints at x — —V%, 
and x = —l/y/2, and a pair of complex conjugate vertical lines bounded by the endpoints 
x ^ -l-v/2±1.0514178378i, and x w -ly/2± 0.38166388451 We have found and additional 
pair of endpoints at x ~ 0.5890850526 ±0.45193582551 There are five pairs of T points; two 
of them are located on the line Rex = — 1\/2. These are x ~ -ly/2± 0.43360353011 and 
x ^ -ly/2± 0.73942467161 The other three pairs are x w -1.0712333535 ± 0.7555078808 i, 
x ^ -1.6123945698±0.80429423591 and x w -0.3186094544±0.9388965869 1 We find four 
bulb-like regions around the latter two pairs of T points. Finally, there is a multiple point 
at x w —0.9415556904, and a complex conjugate pair of multiple points at q — —e ±t7T ^. 
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We have compared the above-described limiting curves with those of a triangular-lattice 
model with free boundary conditions [40]. The agreement is perfect on the whole complex 
x-plane for L = 2, 3, 4. 

9.2 Three-state Potts model (p = 6) 

We have studied the limiting curves given by the sectors Xi,i an d Xi,5> cf- Eq. (7.12). The 
results are displayed in Figure 16. As for the square-lattice case discussed in Section 8.2, 
the limiting curves coincide with those obtained with free boundary conditions in a domain 
containing the real positive f-axis. In particular, the agreement is perfect in the first regime 
Ret> > 0. In the second regime — 1 < Rev < 0, the coincidence holds except on a small region 
close to Ret> = —1, and | Imt> | small for L = 3,4. In both cases, the branches that emerge 
from x = —l/y/3 and penetrate inside the second regime (and defining a closed region), 
change their shape for free boundary conditions (and in particular, the aforementioned closed 
regions are no longer closed). Finally, in the third regime Rev < —1, the limiting curves for 
both types of boundary conditions clearly differ. As for the square-lattice three-state model, 
free boundary conditions usually imply less and smaller closed regions. 

10 Conclusion and outlook 

We have studied the complex-temperature phase diagram of the Q-state Potts model 
on the square and triangular lattices with Q = 4cos 2 (n/p) and p integer. The bound- 
ary conditions were taken to be cyclic so as to make contact with the theory of quantum 
groups [6,7,24,27,28], which provides a framework for explaining how a large amount of the 
eigenvalues of the cluster model transfer matrix — defined for generic values of p — actually 
do not contribute to the partition function Z for p integer. Moreover, for p integer, the exact 
equivalence (2.5) between the Potts and the A p _i RSOS model provides an efficient way of 
computing exactly those eigenvalues that do contribute to Z. Using the Beraha-Kahane- 
Weiss theorem [30] , this permitted us to compute the curves Bl along which partition func- 
tion zeros for cyclic strips of finite width L accumulate when the length N — > oo. 

The RSOS model has the further advantage of associating a quantum number j with each 
eigenvalue, which is related to the number of clusters of non-trivial topology with respect to 
the periodic direction of the lattice and to the spin S z of the associated six-vertex model. 
This number then characterizes each of the phases (enclosed regions) defined by Bl- 

The curves Bl turn out to exhibit a remarkable regularity in L — at least in some respects — 
thus enabling us to make a number of conjectures about the thermodynamic limit L — > oo. 
On the other hand, even a casual glance at the many figures included in this paper should 
convince the reader that the L — > oo limit of the models at hand might well conceal many 
complicated features and exotic phase transitions. Despite of these complications, we ven- 
ture to summarize our essential findings, by regrouping them in the same way as in the list 
of open issues presented in the Introduction: 

1. The points x-fm{Q) and x-{Q) (and for the square lattice also its dual x + {Q)), that act 
as phase transition points in the generic phase diagram, should play a similar role for 
integer p. This can be verified from the figures in which it is more-or-less obvious that 
the corresponding red solid squares will be traversed, or pinched, by branches of Bl in 
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the L — > oo limit. What is maybe more surprising is that also x-qk{Q) has a similar 
property, despite of the profoundly changed physics inside the BK phase. Indeed, in 
most cases, Xbk(Q) is either exactly on or very close to a traversing branch of It 
remains an open question to characterize exactly the nature of the corresponding phase 
transition. 

2. It follows from Conjecture 4.1.2 that for the square lattice, will contain x = — \fQj2 
for p integer and x = —\fQ for Q integer. For the triangular lattice the corresponding 
Conjecture 6.1.2 involves the points x = —2/^/Q for p integer and x = —1/y/Q for 
Q integer. Thus, both lattices exhibit a phase transition on the chromatic line x = 
— 1/y/Q or its dual, but only for integer Q. It is tempting to speculate that the 
chromatic line and its dual might play symmetric roles upon imposing fully periodic 
boundary conditions, but that remains to be investigated. 

3. We have found that with free cyclic boundary conditions, partition functions zeros 
are dense in a substantial region of the phase diagram, including the region \x\ ^> 1. 
See in particular Conjectures 4.3-4.4 for the square lattice and Conjecture 6.4 for the 
triangular lattice. For the Ising model (Q — 2), the finite-size data is conclusive enough 
to make a precise guess as to the extent of that region, cf. Conjectures 7.1-7.2. We have 
argued (in Section 7.4) and observed explicitly (in Sections 8-9) that this feature is 
completely modified by changing to fixed cyclic boundary conditions. Another example 
of the paramount role of the boundary conditions has been provided with the argument 
of Section 8 that when restricting to the sector xi,i one sees essentially the physics of 
free longitudinal boundary conditions. 

4. It is an interesting exercise to compare the limiting curves found here with the numer- 
ically evaluated effective central charge shown in Figs. 23-25 of Ref . [8] . In particular, 
for p = 5 it does not seem far-fetched that the two new phase transitions identified 
in Fig. 23 of that paper might be located exactly at x — —1/y/Q ~ —0.618 and 
x = —y/Q/2 ~ —0.809. These points (for the former point, actually its dual, but we 
remind that the transverse boundary conditions in Ref. [8] are periodic) are among the 
special points discussed in item 2 above. 

5. We have provided some evidence that on the triangular lattice for Q = 4 (i.e., p = oo) 
phases with arbitrary high j will exist close to the point x = — 1. For the square lattice 
we have only found phases with j < 5. This should be compared with the arbitrarily 
high values of S z taken when approaching (Q,x) = (4, —1) from within the BK phase 
in the generic case [7, 34] . 

It would be interesting to extend the study to fully periodic (toroidal) boundary condi- 
tions. This would presumably diminish the importance of finite-size corrections, but note 
that the possibility of the non-trivial clusters having a more complicated topology makes the 
link to the quantum group more subtle. 

Another line of investigation would be to study the Potts model for a generic value of Q, 
i.e., to transpose what we did for the Xi,2j+i to the Ki^j+i- Indeed, studies for v given in the 
complex Q-plane have already been made, for example in Ref. [34] for v = —1, but to our 
knowledge no study exists for Q given in the complex w-plane. Note that the results are very 
different in these two cases. For example, with L fixed and finite, the Beraha number Q = B p 
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are limiting points in the complex Q-plane for fixed v — — 1 (and presumably everywhere 
in the Berker-Kadanoff phase), but v — — 1 is not a limiting point in the complex v- plane 
for fixed Q = B p (p > 4). This is just one example that different limits may not commute 
and the very concept of "a thermodynamic limit" for antiferromagnetic models has to be 
manipulated with great care. 
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A Dimension of the transfer matrix 

The dimension of the transfer matrices J k (p,L) can be obtained in closed form. First 
note that for given p, k — 2j + 1, and L, the dimension of the transfer matrix T k (p, L) 

d[ p \L) = dimT fe (p,L) (A.l) 

equals the number of random walks (with up and down steps) of length 2L steps that start 
at height 1 and end at height k. This random walks have to evolve between a "ceiling" 1 
and a "roof" m — p — 1. 

Let us now proceed in steps. For k — 1 and m = oo we have just the Catalan numbers. 
Thus, if z is the fugacity of a single step, then the ordinary generating function (o.g.f.) is 

/(*) = ''^^ = 1 + £,C L *» (A.2) 

We now keep k — 1, and we introduce the "roof" m. A walk is either empty or consists 
of two independent parts. The first part is between the very first step (necessarily up) and 
the first down step that hits the ceiling (i.e., 1); the second part is the rest of the walk 
(which may be empty). For instance, if p — 4 (m = 3) and L = 3, a possible walk can be 
1-2-3-2-1-2-1. The first part of this walk is 1-2-3-2-1; while the second part of the walk is 
1-2-1. If we take away the first and last steps of the first part (i.e., we are left with 2-3-2), 
we have a walk with m — > m — 1 (as this is equivalent to 1-2-1). Thus, the o.g.f. f(m,z) 
satisfies the equation 

f(m,z) = l + z 2 f(m-l,z)f(m,z) (A.3) 
which is solved by the recurrence 

f(m,z) = — ^ (A.4a) 

1 — z z j(m — 1, z) 

f(l,z) = 1 (A.4b) 

Finally, let us consider the general case with k > 1. In this case, the walk cannot be 
empty, and the first step is necessarily up. There are two classes of walks. In the first one, 



39 



the walk never hits the ceiling 1 again. For instance if p — 4, L — 3, and k = 3, a walk 
belonging to this class is given by 1-2-3-2-3-2-3. So it consists in one step and a walk with 
a raised ceiling (i.e., 2-3-2-3-2-3 is equivalent to 1-2-1-2-1-2 with roof m = 2). In the 
second class, the walk does hit the ceiling somewhere for the first time, so we can split the 
walk into two independent parts as in the preceding paragraph. Thus, the o.g.f. satisfies the 
equation 

f(m,k,z) = z f(m — 1, k — 1, z) + z 2 f(m — 1, z) f(m, k, z) (A. 5) 

which can be solved by the recurrence 

f( i \ zf(m-l,k- l,z) 

f(m,k,z) = — — - A.6a 

1 — z 2 j(m — 1, z) 

/(m,M) = f(m,z) (A.6b) 
where f(m,z) is given by (A. 4). The dimensions d^\L) can be read off immediately 

oo 

f(m,k,z) = Y,d { k\L)z 2L (A.7) 

L=0 

In the particular case p = 4, we easily find that 

/(3,M) = = l + ^-V^ (A.8a) 

.2 



L=l 



oo 

/(3,3,,) = = ^2 L -^ 2i (A.8b) 

L=l 

For the other cases, we can get closed formulas for the generating functions, and obtain the 
result 

4 P) (L) = hn P+j (L) - 7 (B+ i)p-w(L)) • (A.9) 

n>0 

where k = 2j + 1 and we have defined Jj(L) = for j > L. The 7j(L) are given by 

^-a")-a"-o-^i(w)- < aio) 

This result can also be obtained by another method, which consists of calculating the number 
7j(L) of states of highest weight with spin S = S z = j for the vertex model and taking into 
account the coupling of U q (SU (2)) between different j for p integer [28]. Yet another method 
consists in relating to the number of paths on the Dynkin diagram going from 1 to 
2j + 1 and using the eigenvectors of the adjacency matrix [27]. 

The 7j(L) can also be interpreted as the dimension of the transfer matrix in the FK 
representation with j bridges, i.e., for a generic (irrational) value of p. In that context, 
Eq. (A.9) represents the reduction of the dimension that takes case at p integer when going 
from the FK to the RSOS representation (with spin j), and thus, is completely analogous 
to Eq. (2.6) for the generation functions. 

On the chromatic line x = —1/y/Q, jj(L) is replaced by a smaller dimension Tj(L), 
because the operator V — Y[ Vi is a projector (V 2 = V) that projects out nearest-neighbor 
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connectivities (i.e. the action of V on states with nearest neighbours connected gives zero). 
We do not know of any explicit expression for Tj(L), but it verifies the following recursion 
relation [47] 



with the convention that ^j(L) = for j < and the conditions r?(L) = for j > L, 
T L (L) = 1, and r (l) = 0. In particular, it can be shown that T (L) = M L _ 1; where 
M L _i is a Motzkin number and corresponds to the number of non-crossing non-nearest 
neighbor partitions of {1, . . . , L} (i.e., it is the dimension of the cluster transfer matrix in 
the case of free longitudinal boundary conditions and x = —1/yfQ). Note that in the RSOS 
representation (in the case of p integer), we cannot reduce the dimension of the since 
although V is a projector in the RSOS representation too the states which are projected 
out are linear combinations of the basis states (corresponding to a given configuration of the 
heights), and not simply basis states as in the case of the FK representation. But because 
of Eq. (2.6), the number dj(L) of non null eigenvalues of T 2 j +1 is given by Eq. (A. 9) with 
7j(L) replaced by Tj(L). In particular, for Q = 3 we find using the recursion relation that 
di(L) = d 5 (L) = 2 L ~ 2 and d 3 (L) = 2 L_1 . Indeed, for x = -l/^/Q, the three-state Potts 
model is equivalent to a homogeneous six-vertex model with all the weights equal to 1 [19] 
(note that this six- vertex model is different from the one we considered before). 



r (£ + i) 
r,-(L + i) 



Fj-i(L) + Tj(L) + T j+1 (L) for j > 



(A.lla) 
(A.llb) 
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Table 1: First L coefficients for the leading eigenvalue A^i(L) coming from the sector xi,i 
for a square-lattice strip of width L. 
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Table 2: First min(2L — 1, 7) coefficients b k for the leading eigenvalue A^i(L) coming from 
the sector for a triangular-lattice strip of width L. 
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Square— lattice Potts model 




Q 

Figure 1: Generic phase diagram for the two-dimensional Potts model in the (Q, t>)-plane. 
The solid black curve in the ferromagnetic (v > 0) region shows the standard ferromagnetic 
phase transition curve vfm(Q), and the blue dashed curve is its analytic continuation vbk(Q) 
into the antiferromagnetic region. This latter curve acts as an RG attractor for the Berker- 
Kadanoff phase (the orange hatched region). This is separated from the limit of infinite 
temperature (red dashed curve) by the antiferromagnetic phase-transition curve v + (Q) (solid 
black curve in the v < region), and from the v — > — oo limit by its counterpart V-(Q) 
(dot-dashed blue curve). The red horizontal dotted curve represents the zero-temperature 
antiferromagnet (v = —1). The pink vertical lines show the Beraha numbers Q = 4 cos 2 (7r/p) 
(p = 2, 3, . . .): the phase diagram on these lines is different from the generic one shown here 
and forms the object of the present article. Note that the exact functional forms of the curves 
vfm(Q), vbk(Q), and v±(Q) are lattice-dependent; the figure shows their explicit forms for 
the square-lattice model. 
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_ 1 _ 2 l/2 _ 2 l/2 _j _ 2 -l/2 j.21/2 q j 

Figure 2: Phase diagram and RG flows for the Q 2 state model (top) and the Q = 2 Ising 
model (bottom), on the real x-axis. Filled (resp. empty) circles correspond to critical (resp. 
non-critical) fixed points. 




Figure 3: RSOS lattice (solid thick lines) and label convention for the basis in the height 
space for a square-lattice of width L = 2 (dashed thinner lines). The thick black arrow shows 
the transfer direction (to the right). 
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Figure 4: RSOS lattice (solid thick lines) and label convention for the basis in the height 
space for a triangular- lattice of width L = 2 (dashed thinner lines). The thick black arrow 
shows the transfer direction (to the right). 
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Square lattice p=4 L=2 F Square lattice p=4 L=3 F 




Figure 5: Limiting curves for the square-lattice RSOS model with p = 4 and several widths: 
L — 2 (a), L — 3 (b), and L = 4 (c). For each width L, we also show the partition- 
function zeros for finite strips of dimensions Lp x (10L)p (black □), Lp x (20L)p (red o), and 
Lp x (30L)p (brown A). Figure (d) shows all these limiting curves together: L = 2 (black), 
L = 3 (red), L = 4 (green). The solid squares ■ show the values where Baxter found the 
free energy. The symbol x in (a) marks the position of the found isolated limiting point. In 
the regions displayed in gray (resp. white) the dominant eigenvalue comes from the sector 
Xi,3 (resp. The dark gray circles correspond to (1.5) 
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Square lattice p = 5 L=2 F Square lattice p = 5 L=3 F 




Figure 6: Limiting curves for the square-lattice RSOS model with p = 5 and several widths: 
L — 2 (a), L — 3 (b), and L — 4 (c). Figure (d) shows all these curves together: L = 2 
(black), L = 3 (red), L = 4 (green). The solid squares ■ show the values where Baxter found 
the free energy. In the regions displayed in light gray (resp. white) the dominant eigenvalue 
comes from the sector xi,3 (resp. The dark gray circles correspond to (1.5) 
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Square lattice p=6 L=2 F 



Square lattice p=6 L=3 F 




Square lattice p=6 L=4 F 



(b) 

Square lattice p=6 




Figure 7: Limiting curves for the square-lattice RSOS model with p = 6 and several widths: 
L = 2 (a), L = 3 (b), and L = 4 (c). Figure (d) shows all these curves together: L = 2 
(black), L = 3 (red), L = 4 (green). The solid squares ■ show the values where Baxter 
found the free energy. In the regions displayed in light gray (resp. white) the dominant 
eigenvalue comes from the sector xi,3 (resp. 111 the regions displayed in a darker gray 

the dominant eigenvalue comes from the sector Xi5- The dark gray circles correspond to 
(1.5) 
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Square lattice p=°° L=2 F Square lattice p=°° L=3 F 




Figure 8: Limiting curves for the square-lattice RSOS model with p = oo (Q = 4) and 
several widths: L — 2 (a), L — 3 (b), and L — 4 (c). Figure (d) shows all these curves 
together: L = 2 (black), L = 3 (red), L = 4 (green). The solid squares ■ show the values 
where Baxter found the free energy. In the regions displayed in light gray (resp. white) the 
dominant eigenvalue comes from the sector xi,3 (resp. Xi,i)- m ^ ne re gi° ns displayed in a 
darker gray the dominant eigenvalue comes from the sector xi,5- 
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Triangular lattice p=4 L=2 F Triangular lattice p=4 L=3 F 




Figure 9: Limiting curves for the triangular-lattice RSOS model with p = 4 and several 
widths: L = 2 (a), L = 3 (b), and L = 4 (c). For each width L, we also show the partition- 
function zeros for finite strips of dimensions L F x (10L)p (black □), L F x (20L) P (red o), and 
Lp x (30L)p (brown A). Figure (d) shows all these limiting curves together: L = 2 (black), 
L = 3 (red), L = 4 (green). The solid squares ■ show the values where Baxter found the 
free energy. The symbol x in (a) marks the position of the found isolated limiting point. In 
the regions displayed in gray (resp. white) the dominant eigenvalue comes from the sector 
Xi,3 (resp. The gray ellipse corresponds to (Rex + l/\/2) 2 + 3(Imx) 2 = 3/2. This 

curve goes through the points x = —e ±t7T ^. 
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Triangular lattice p=5 L=2 F 



Triangular lattice p=5 L=3 F 
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Figure 10: Limiting curves for the RSOS model with p = 5 and several widths: L = 2 (a), 
L = 3 (b), and L = 4 (c). Figure (d) shows all these curves together: L = 2 (black), L = 3 
(red), L = 4 (green). The solid squares ■ show the values where Baxter found the free 
energy. In the regions displayed in light gray (resp. white) the dominant eigenvalue comes 
from the sector Xi i3 (resp. 
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Triangular lattice p=6 L=2 F Triangular lattice p=6 L=3 F 




Figure 11: Limiting curves for the triangular-lattice RSOS model with p = 6 and several 
widths: L — 2 (a), L — 3 (b), and L — 4 (c). Figure (d) shows all these curves together: 
L = 2 (black) , L = 3 (red) , L = 4 (green) . The solid squares ■ show the values where Baxter 
found the free energy. In the regions displayed in light gray (resp. white) the dominant 
eigenvalue comes from the sector xi,3 (resp. 111 the regions displayed in a darker gray 

the dominant eigenvalue comes from the sector Xi,5- 
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Triangular lattice p=°° L=2 F Triangular lattice p=°° L=3 F 




Figure 12: Limiting curves for the triangular-lattice RSOS model with p = oo (Q = 4) and 
several widths: L — 2 (a), L — 3 (b), and L = 4 (c). Figure (d) shows all these curves 
together: L = 2 (black), L = 3 (red), L = 4 (green). The solid squares ■ show the values 
where Baxter found the free energy. In the regions displayed in light gray (resp. white) 
the dominant eigenvalue comes from the sector xi,3 (resp. m regions displayed in a 

darker gray the dominant eigenvalue comes from the sector xi ,5- in (c), an even darker gray 
marks the regions with a dominant eigenvalue coming from the sector Xij- The gray circle 
corresponds to (7.8). 
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Square lattice p=4 L=2 X 



Square lattice p=4 L=3 X 




Square lattice p=4 L=4 X 



Square lattice p=4 




Figure 13: Limiting curves for the square-lattice RSOS model with p = 4 and several widths: 
L — 2 (a), L — 3 (b), and L — 4(c) when only the sector xi,i is taken into account. Figure (d) 
shows all these curves together: L = 2 (black), L = 3 (red), L = 4 (green). The solid squares 
■ show the values where Baxter found the free energy. The dark gray circles correspond to 
(1.5) 
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Square lattice p=6 L=2 X 



Square lattice p=6 L=3 X 




Square lattice p=6 L=4 X 



Square lattice p=6 




Figure 14: Limiting curves for the square-lattice RSOS model with p = 6 and several widths: 
L — 2 (a), L — 3 (b), and L — 4(c) when only the sectors xi,i an d Xi,5 are taken into account. 
In the regions displayed in dark gray (resp. white) the dominant eigenvalue comes from the 
sector xi,5 (resp. Figure (d) shows all these curves together: L = 2 (black), L = 3 

(red), L = 4 (green). The solid squares ■ show the values where Baxter found the free 
energy. The dark gray circles correspond to (1.5) 
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Triangular lattice p=4 L=2,j 



Triangular lattice p=4 L=3,j 




Triangular lattice p=4 L=4,, 



Triangular lattice p=4 




Figure 15: Limiting curves for the triangular-lattice RSOS model with p = 4 and several 
widths: L — 2 (a), L — 3 (b), and L — 4 (c) when only the sector is taken into account. 
Figure (d) shows all these curves together: L = 2 (black), L = 3 (red), L = 4 (green). The 
gray ellipse corresponds to (Rex + l/v^) 2 + 3(Imx) 2 = 3/2. This curve goes through the 
points x = —e ±ln / 4 . 
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Triangular lattice p=6 L= 2 X Triangular lattice p=6 L— 3 X 
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Figure 16: Limiting curves for the triangular-lattice RSOS model with p = 6 and several 
widths: L — 2 (a), L — 3 (b), and L = 4 (c) when only the sectors an d Xi,5 are taken 
into account. In the regions displayed in dark gray (resp. white) the dominant eigenvalue 
comes from the sector xi,5 (resp. Figure (d) shows all these curves together: L = 2 

(black), L = 3 (red), L = 4 (green). The solid squares ■ show the values where Baxter 
found the free energy. 
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